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1.  INTRODUCTION 


There  are  two  main  types  of  vibratory  forties  induced  by  a propeller. 

The  first  type  is  known  as  bearing  forces  , and  consists  of  induced  vibratory 
forces  on  the  propeller  caused  by  the  incoming  non-uniform  wake  of  the  ship. 
These  forces  are  delivered  to  the  hull  shrough  the  stem  bearing.  The 
second  type,  known  as  surface  forces , consists  of  induced  vibratory  loads 
through  water  on  a structure  near  the  propeller  such  as  the  hull  or  the 
rudder;  They  are  caused  by  time  dependent  variations  of  the  flow  field,  due 
to  the  rotation  of  the  propeller.  As  a result,  these  latter  forces  may 
arise  even'  if  the  propeller  operates  in  a steady  flow. 

One  should  note  at  this  point  that  the  term  forces  is  used  in  this 
context  in  a broader  sense  and  should  be  interpreted  as  meaning  any 
type  of  hydrodynamic  excitation  arising  from  the  propeller  whether  this 
is  a force  or  a moment. 

It  is  clear  from'  the  above  definitions  that  in  the  context  of  this 
report  only  surface  forces  are  of  interest . The  above  subdivision  of 

t 

propeller  forces  is  by  no  means  ■universal.  Different  authors  use  dif- 
ferent terms.  F.M. Lewis  [ 2,3  ],  for  example  , distinguishes  between 

surface  (acting  on  the;  hull  surface  excluding  the  rudder)  and  rudder  forces. 

A further  distinction  is  often  made  among  surface  forces  by  subdividing 
them  into  pitch  unbalance  forces  and  blade  frequency  forces. 

Had  we  had  a perfectly  shaped  propeller  operating  at  constant  r.p.m. 
in  a uniform  axial  Wake  then  symmetry  would  require  the  generation  of 
constant  torque,  constant  thrust , zero  lateral  force  variation,  and-  zero 
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moment  variation  about  any  axis  normal  to  its  shaft  axis. 

Pitch  unbalance  forces  and  moments  are  generated  by  non-uniformities 
embodied  in  the  propeller  blades  during  manufacturing  or  due  to  blade 
damage  occurred  as  a result  of  the  propeller  running  into  an  obstruction. 

The  effect  of  such  irregularities  is  the  introduction  of  first  order 
(frequency  same  as  that  of  the  propeller  shaft)  forces  and  moments  at  the 
propeller  shaft  bearing.  If  the  propeller  operates  in  a non-uniform  wake 
as  is  usually  the  case  with  surface  ships.  Ditch  unbalance  would  also  result 
in  the  introduction  of  superimposed  variation  effect  with  a first  order 
fundamental  frequency  but  which  may  have  harmonic  components  depending 
on  which  frequencies  are  excited  by  the  non-uniform  flow. 

Blade  frequency  forces  exist  regardless  of  the  degree  of  perfection 
in  the  manufacture  of  the  propeller.  As  one  would  expect  they  vary  at 
blade  frequency  (equal  to  propeller  r.p.m.  x no.  of  blades)  or  at  one  of 
its  harmonics. 

The  hull  surfaces  on  which  surface  forces  act  are  classified  as  tip- 
clearance  s'M'faces  and  axial  clearance  surfaces.  In  single  screw  ships  tip 
cl Barance  surfaces  are  the  hull  surfaces  above  the  propeller  whilst  axial 
clearance  surfaces  include  struts,  bossings,  and  the  rudder. 

The  operation  of  the  propeller  in  a non-uniform  wake  tends  to  increase 
the  surface  forces  above  the  level  they  would  have,  had  the  inflow  been 
uniform.  These  wake-induced  force  movements  are  known  as  wake -reflection 
forces . 

The  rudder  in  a single  screw  ship  extends  inside  the  vortex  sheet 
flow  behind  the  propeller,  the  various  components  of  which  produce  vibratory 
forces  of  frequency  equal  to  propeller  blade  frequency. 
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The  study  of  propeller-induced  rudder  vibration  has  much  more  recently 
been  developed.  Its  importance  has  been  augmented  as  ship  speeds  and 
propeller  thrust  loadings  have  increased.  Most  of  the  relevant  research 
has  been  conducted  over  the  last  decade  and  has  frequently  been  related 
to  the  field  of  hydroelasticity. 

Vibration  of  the  rudder  above  a certain  level  may  either  cause  local 
failure  or  may  in  turn  excite  vibration  of  the  entire  hull  ("spring  mass 
effect").  Prediction  of  dynamic  response  characteristics  of  the  rudder 
to  hydrodynamic  loads  can  be  of  great  assistance  to  the  designer. 

Although  the  present  experiments  carried  out  at  M.I.T.'s  water 
tunnel  were  trying  to  simulate  conditions  behind  a particular  ship,  the 
results  may  be  valid  for  other  vessels  too. 

The  ship  in  question,'  the  Sea  Control  Ship,,  is  a single 
screw,  high  speed,  transom  stem  type  ship.  It  also  comprises 
the  additional  feature  (which  is  getting  quite  common  nowadays)  of 
having  the  design  location  of  the  rudder  moved  slightly  to  the  starboard 
of  the  center  line  of  the  propeller  shaft  axis . 

The  reason  behind  this  is  twofold.  It  causes  the  rudder  to  miss 
the  hub  vortex  of  the  propeller  and  enables  the  propeller  shaft  to  be 
removed  without  heeding  to  unship  the  rudder  every  time  this 
is  desired. 

In  such  a case  maintenance  or  repair  time  can  be  reduced  considerably 
which  is  an  important  advantage  to  ships  with  high  availability  requirements . 
As  a result,  one  of  the  objects  of  the  experiment  was  to  examine  whether  this 


4 


special  feature  has  a detrimental  effect  on  the  level  of  induced  vibration 
or  not. 

One  may  note,  however,  that  this  feature  which  is  getting  quite  usual 

practice  on  single  screw  ships  is  even  more  common  to  twin  screw  ships 

where  the  two  rudders  are  not  aligned  with  the  propeller  axes.  Thus,  the 

results  of  this  experiment  can  be  of  some  use  to  twin  screw  ships  by  giving 

some  indication  of  the  effect  on  vibrations  of  various  transverse  positions. 

* 

Since  longitudinal  propeller-rudder  clearance  is  a critical  parameter 
in  determining  vibratory  force  levels  on  the  rudder,  it  is  important  to 
establish  a precise  definition  at  the  outset.  In  accordance  with  current 
NAVSEC  practice,  this  clearance  is  defined  as  the  distance  between  the  center 
of  the  blade  when  it  is  at  its  uppermost  position  and  the  leading  edge 
of  the  rudder  measured  on  a lice  parallel  to  the  propeller  shaft  at  a 
distance  equal  to  70%  of  the  propeller  radius.  This  definition  differs 
from  that  used  by  F.  M.  Lewis  [2],  [3],  where  the  distance  is  measured  at 
the  shaft  centerline,  rather  than  at  the  70%  radius. 
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2.  SKj.?  SIMULATION  IN  THE  WATER  TUNNEL 

Since  in  our  experiment  we  were  trying  to  investigate  what  happens 
behind  a given  ship,  it  was  necessary  for  us  to  make  maximum  use  of  the 
information  we  had  been  given  on  the  ship's  design  features  and  its  opera- 
tional characteristics. 

To  facilitate  our  task,  the  project  sponsor  provided  us  with  drawings 
of  the  stern  part  of  the  ship,  specifying  the  limits  of  the  permissible 
changes  of  rudder  longitudinal  and  transverse  positions.  The  ship  operating 
conditions  were  also  given  to  us  in  non-dimensional  form  (by  providing 
values  of  J,  a,  and  Kj.)  together  with  the  results  of  the  ship's  wake 
survey . 

In  addition  to  all  this  information,  we  were  given  the  model  of  the 
propeller.  This  immediately  determined  the  model/ship  scale  ratio,  ap- 
proximately .0410,  and  took  care  of  the  requirement. 

To  be  able  to  satisfy  both  the  J and  o requirements  we  made  use  of 
dimensional  analysis  (see  Appendix  V)  and  came  forward  with  a set  of 
tunnel  operating  conditions'  options.  Deciding  which  of  these  options  was 
to  be  selected  was  left  to  the  results  of  the  natural  frequency  analysis. 

Combining  all  the  above  information  as  well  as  information  obtained 
from  close  examination  of  the  ship's  drawings,  the  design  requirements  of 
the  model  may  be  summarized  as  follows: 

1)  The  model  must  satisfy  the  rudder  lateral  movement  requirements.  In 
other  words,  the  rig  must  be  such  as  to  enable  the  accurate  location  of 
the  dynamometer-rudder  subassembly  at  exactly  the  seven  predetermined 
positions. 


6 


2)  It  must  also  be  able  to  simulate  the  slight  flow  inclination  (~4.5°) 
relative  to  the  "vertical"  (i.e.  normal  to  baseline)  rudder  stock  (here 
modelled  by  the  dynamometer) . At  present,  however,  the  tunnel's  propeller 
drive  system  does  not  allow  for  changes  in  the  shaft  angle.  As  a result, 
the  propeller  shaft  must  be  horizontal  and  the  rudder  dynamometer  must  bear 
the  inclination.  Either  way  would  produce  similar  results. 

3)  The  supporting  rig  must  permit  changes  in  the  angle  of  attack  (rudder 
deflection) . 

4)  The  model  must  allow  for  variations  in  the  axial  propeller-to-rudder 
clearance.  One  should  note  that  "clearance'''  here  is  defined  as  the  distance 
between  the  centre  of  a blade,  where  it  is  at  its  uppermost  positions,  and 
the  leading  edge  of  the  rudder  measured  on  a line  parallel  to  the  propeller 
shaft  centerline  and  at  a distance  equal  to  70%  of  the  propeller  radius,  R. 

5)  It  must  also  be  made  watertight  whenever  leak  paths  are  anticipated. 

6)  The  supporting  rig  must  be  rigid  enough  to  allow  for  very  small  axial 
and  lateral  displacement  of  the  assembly. 

7)  The  model  must  incorporate  certain  appendages  (in  the  form  of  screens) 
which  will  simulate  the  ship's  axial  wake  as  closely  as  possible. 

8)  All  positions  and  angular  changes  must  involve  a minimum  number  of 
adjustments. 

9)  Finally,  the  model  must  have  a natural  frequency  lying  above  the  an- 
ticipated operating  blade  frequency. 

A description  of  the  final  version  of  the  model  is  given  ir.  the 
following  section.  The  various  considerations  that  led  to  It  however, 
are  discussed  in  more  detail  in  Appendix  I. 
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Description  of  the  Model 

A view  of  the  model  is  given  in  fig.(l).  The  aluminum 
window  is  bolted  down  to  position  on  the  tunnel.  A thick  steel  plate 
(the  "infinite  mass")  is  bonded  unto  the  window  by  means  of  a 9/16"  thick 
layer  of  ECCOBOND  epoxy  mixture  (2  parts  ECCOBOND  to  3 parts  of  catalyst 
by  weight).  This  method  of  attaching  plate  to  window  was  chosen  as  our 
final  attempt  to  decrease  the  system’s  natural  frequency  by  lowering  its 
stiffness. 

The  dynamometer  is  held  into  position  by  a massive  aluminum  clamp 
bolted  on  top  of  the  steel  plate.  Lateral  movement  can  be  achieved  by 
allowing  the  clamp  to  move  inside  a groove  cut  out  in  the  plate.  The 
clamp  must  be  sufficiently  large  to  cover  the  entire  groove  at  all  positions. 
Note  that  the  groove  is  large  enough  to  ensure  access  to  1/2  the  lateral 
positions  required.  The  other  half  can  be  obtained  by  rotating  the  window 
180°.  This  method, although  time  consuming,  has  two  great  advantages. 

It  wastes  a minimum  amount  of  steel  mass  and  requires  a smaller  clamp. 

The  rudder  is  held  tight  against  the  tapered  bottom  end  of  the  sensor. 

A "tiller"  used  for  angular  adjustments  is  fixed  on  top  of  the  dynamometer. 
Whenever  a change  in  angle  of  attack  is  to  be  made,  the  clamp  is  slightly 
loosened  and  the  rudder  assembly  is  rotated  to  its  desired  position.  A 
tapered  pin  sliding  through  two  holes  on  the  tiller  and  the  clamp,  prevents 
any  further  movement. 

Sealing  is  provided  by  an  0-ring  between  the  steel  plate  and  the 
clamp  and  another  0-ring  and  a gasket  between  the  jaws  of  the  clamp. 


To  further  simulate  the  ship's  hull  by  providing  the  proper  propeller 


tip  clearance,  it  was  found  necessary  t:  use  a splitter  plate.  This  is 


i'if'.l:  Assembly  drawing  of  equipment  mounted  in  the  water  tunnel  test  section 
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held  into  position  1 3/4  inches  from  the  top  of  the  test  section  by  means 
of  an  aluninum  streamlined  boat.  Like  the  steel  plate,  both  splitter  plate 
and  boat  have  a transverse  groove  cut  into  them  to  allow  for  sideways  move- 
ment of  the  clamp-dynamometer  subassembly.  The  splitter  plate  extends 
across  the  whole  width  of  the  test  section  leaving  practically  zero-clear- 
ance at  the  sides  to  ensure  the  complete  splitting  of  the  flow  above  and 
under  the  plate. 

Although  the  splitter  plate  is  longer  than  the  previous  one  used  for 
other  experiments  in  the  tunnel[9],  it  still  extends  less  upstream  than 
the  test  section  wall  and  thus  manages  to  reduce  the  boundary  layer  thick- 
ness (see  Appendix  III),  finally,  to  avoid  flow  separation,  the  leading 
edge  was  faired  as  shown  in  fig.(l). 

The  groove  in  the  splitter  plate  is  covered  by  means  of  a set  of  inter- 
changeable coverplates,  one  of  which  carries  the  rudder  skeg  (fixed  part 
of  the  rudder).  The  skeg  proportions  had  to  be  altered  slightly  to  make 
it  possible  for  the  dynamometer  to  go  through  it  with  clearance- 

The  vee-strut  assembly  slipping  over  the  propeller  shaft  is  used  to 
simulate  the  strut  on  the  ship.  Three  set  screws  holding  it  tightly  against 
the  drive  shaft  housing  make  sure  that  it  remains  at  a constant  distance 
from  the  propeller  as  the  longitudinal  position  of  the  latter  is  changed. 

finally,  to  simulate  the  non-uniform  wake  behind  the  ship,  it  was  de- 
cided to  introduce  in  the  flow  a wake  screen,  attached  near  the  leading 
edge  of  the  splitter  plate,  capable  of  creating  a velocity  profile  similar 
to  that  of  the  ship.  The  variable  solidity  of  the  screen  as  n function  of 
distance  from  the  propeller  shaft  centreline  was  developed  in  a wake  survey 
carried  out  in  the  tunnel  described  in  Appendix  HI. 
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3 . EXPERIMENTAL  PROCEDURE 

Although  we  were  primarily  interested  in  examining  the  effect  of  pro- 
peller to  rudder  orientation  on  the  vibratory  loading  of  the  rudder  (dynamic 
tests),  it  was  found  worthwhile  to  extend  the  investigation  a step  further 
to  include  the  effect  on  the  steady  state  characteristics  of  both  propeller 
and  rudder  (steady-state  tests). 

This  decision  was  taken  in  order  to  ensure  that  the  longitudinal  and 
lateral  position  of  the  rudder,  with  respect  to  the  propeller,  the  dynamic 
tests  would  indicate  as  being  the  "optimal",  would  not  have  a detrimental 
effect  on  the  steady  state  performance  characteristics. 

Thus,  the  overall  experiment  was  divided  into  two  main  parts,  the  first 
dealing  with  dynamic  tests  and  the  second  with  steady  state  tests.  For  sim- 
plicity these  are  treated  here  separately. 

Dynamic  Tests  - Modifications  to  the  Original  Plan 

By  the  time  we  started  the  first  tests , in  early  July  1974 , we  had  already 
decided  to  modify  our  original  testing  plans.  It  was  then  that  we  decided 
to  examine  only  the  centreline  and  starboard  lateral  positions  of  the 
rudder, based  on  indications  obtained  after  the  first  test  runs,  that  the 
optimum  lateral  position  was  near  the  design  value  (which  corresponds  with 
1"  to  the  starboard  on  the  model).  This  decision  was  also  enforced  by 
another  consideration.  Obtaining  the  port  lateral  positions  (as  already 
described)  meant  rotating  the  window  by  180°  and  refastening  it 
on  the  boat.  Because  of  the  absence  of  symmetry  in  the  shape  of  the 
boat  this  last  step  would  inevitably  involve  drilling  a new  set  of 
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holes  into  the  window,  thus  increasing  ths  nurahar  of  possible  air-leak 
paths  into  the  system  which  we  knew  woul*.*'  hinder  our  experiment. 

After  a number  of  trial  runs,  we  discovered  that  there  was  little 
advantage  gained  by  using  a smaller  increment  ,*f  longitudinal  movement 
of  the  propeller  near  the  design  position  (at  a clearance  ratio  of  about  .57). 
Instead,  one  centimeter  intervals  were  used  throughout  the  travel  of  the 
propeller  with  the  only  exception  that  in  the  three  last  points  (at  very 
large  clearance  ratios)  the  increment  was  increased  to  r.:o  cms. 

The  final  change  involved  the  investigation  of  the  effect  of  changes 
in  angle  of  attack  on  the  vibratory  loading  of  the  rudder.  For  reasons 
explained  in  Appendix  II,  we  were  not  able  to  get  any  results  at  angles 
of  attack  greater  than  5°.  As  a result,  our  original  plans  for  a thorough 

investigation  of  this  effect  had  to  be  abandoned. 

Dynamic  Test  Procedure 

With  the  rudder,  at  one  of  the  four  lateral  locations  to  be  investi- 
gated, the  tunnel  was  first  run  at  a lew  water  speed  and  a static  pressure 
close  to  the  desired  testing  value  (236-nmHg)  in  an  effort  to  drive  as 
much  air  out  of  solution  as  possible.  The  air  thus  released  was  gathered 
in  the  vacuum  chamber  at  the  top  of  the  tunnel,  and  upon  further  reductions 
•i  pressure,  was  driven  completely  out  of  the  system. 

After  a period  of  time  varying  between  20  minutes  and  1 hour,  both  the 
propeller  r.p.m.  and  the  tunnel  water  speed  were  increased  until  the  desired 
testing  values  were  attained  (1500  r.p.m.  and  17.3  ft/sec  respectively). 

Some  additional  time  was  then  allowed  for  the  tunnel  to  reach  equilibrium 
at  these  new  conditions  during  which  frequent  adjustments  were  necessary. 
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Once  equilibrium  was  reached,  a fact  which  was  verified  by  the  lack  of 
fluctuations  in  the  level  of  the  "blue  fluid"  inside  the  manometer,  we 
were  able  to  start  taking  readings. 

For  all  the  dynamic  tests  involving  0°  angle  of  attach  the  dynamometer 
amplifiers'  gain  was  set  at  its  ntsdimjm  value.  This  was  done  in  ant ieipat ion  of 
the  fact  that  the  expected  steady  state  loads  at  0°  rudder  angle , were  not  large 
enough  to  cause  saturation  of  the  amplifiers.  Nevertheless,  the  two  dynamaneter 
response  signals  F anc  B„,  were  examined  prior  to  being  averaged,  were  viewed 
on  an  oscilloscope  (TEKTRONIX  model  561B)  and  a biasing  voltage  was  applied 
whenever  necessary.  After  this  check  for  saturation,  the  two  response  signals 
were  supplied,  one  at  a time  (via  a six-way  switch),  to  the  waveform  eductor 
"signal  input"  jack.  A detailed  description  of  the  dynamometer  and  the 
waveform  eductor  is  given  in  Appendix  II.) 

A switch  on  the  eductor’s  front  panel  is  used  to  select  one  of  the 
three  possible  operating  modes  of  the  instruments.  In  tne  "analyze  and 
readout"  mode  the  instrument  simultaneously  samples  the  input  signal,  stores 
the  repetitive  synchronized  component  of  the  signal  in  its  memory  and  sup- 
plies the  contents  thereof  to  the  output.  With  a characteristic  time  con- 
stant, T,  set  at  20  secs,  it  took  the  eductor  approximately  2 minutes  to 
analyze  and  store  the  input  signal.  At  the  end  of  the  2 minutes,  the  mooe 
switch  was  turned  to  the  "readout  only"  position.  In  this  mode  no  further 
sampling  of  the  '.r.oct  sL~-  a.:  lakes  place,  but  the  instrument  simply  scans 
its  mericry  and  supplies  . re.  contents  thereof  to  the  output.  Sucn  a signal 
was  uc"d  to  drive  the  oscilloscope  which  was  triggered  externally  with  tne 
"gate  signal",  wh_.  ■ i,  a positive  pulse  available  for  the  duration  of  each 
sweep,  at  the  gate  jaci1  c-n  the  eductor's  front  panel.  An  oscilloscope 
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camera  (TEKTRONIX  model  C-59)  was  then  used  to  make  a permanent  record 
of  the  output  response.  After  the  picture  was  taken,  the  signal  was 
erased  (by  turning  the  mode  switch  to  the"erase"  position)  and  then  the 
whole  cycle  was  repeated  for  the  second  dynamometer  signal. 

In  order  for  the  eductor  averaging  sweeps  to  be  synchronized  with  the 
rotation  of  the  propeller  an  external  trigger  pulse  had  to  be  supplied  to 
the  "trigger  in"  jack.  This  was  accomplished  using  a magnetic  pick-up, 
located  at  the  rear  end  of  the  propeller  drive,  which  activated  a - ke 
generator  every time  a tooth  of  the  five- tooth  gear  rotati**  • shaft, 

parsed  over  It.  With  a five-tooth  gear  and  an  appropriate  sv.  juration 

(about  7.8msecs)  it  can  be  arranged  so  that  the  eductor  examines  the  input 
signal  five  times  per  propeller  revolution.  This  provides  a better  system 
response  criterion  since  the  eductor  will,  at  each  sweep,  integrate  the 
input  signal  with  a different  blade  at  a particular  position  so  that  indi- 
vidual blade  non-uniformities  will  be  removed  in  the  averaging  process. 

One  should  note  that  the  two  pulses,  the  one  triggering  the  eductor  and 
the  one  available  at  the  "gate"  jack  are  perfectly  synchronised. 

After  one  full  run  was  completed  (there  were  approximately  IT  longi- 
tudinal positions  per  run)  a check  for  repeatability  was  undertaken.  This 
was  in  fact  nothing  more  than  a repetition  of  the  readings  previously  taken, 
at  randonly  chosen  points.  It  was  during  one  of  these  checks  that  we  dis- 
covered how  amplifier  saturation  could  affect  our  readings  as  well  as  the 
Importance  of  high  air  content  . s the  water.  After  we  satisfied  ourselves 
that  the  readings  originally  taken  could  be  reproduced  we  would  go  on  and 


change  the  lateral  position  of  the  rudder. 


Finally,  *s  startK  testily  for  dasss  in  sflc  o?  attack  (naASer 
wflectiai  angle).  Tot*  these  tests,  d»  to  the  increase  in  the  Iced,  the 
amplifier  gain  had  to  be  GrcmaC  re  Its  tJrisLss  wait*.  Nevertheless,  jatura- 
ticnof  the  ssUfir.s  saie  if  Isgossible  for  'is  to  get  results  for  nsd&r  aigles 
greater  than  S®.  ths  checked  the  effttt  on  the  vibratory  loading  i changing 
the  angle  frc *a  -5®  to  *S®  at  no  cifferem  1 cegitudinal  cl«ra^  at  il» 
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In  the  first  part , the  program  makes  use  of  the  results  obtained  in 
the  dynamic  calibration  (see  Appendix  IV)  to  develop  a matrix  of  the  com- 
plex influence  coefficients.  This  matrix  is  then  used  in  the  second  part 
of  the  program  in  the  reduction  of  the  data  obtained  during  the  experiment. 
The  final  output  is  given  in  non-dimensional  form  as  force  and  moment 
coefficients.  In  both  parts  of  the  program  the  input  consists  of  the ‘ampli- 
tude and  phase  values  of  the  recorded  signals  of  the  two  channels  obtained 
at  a given  point. 

The  amplitude  of  the  signal  is  equal  to  half  the  peak  to  peak  value 
directly  measured  on  the  photograph  multiplied  by  the  sensitivity  of  the 
scale. 

The  phase,  on  the  other  hand,  is  taken  as  the  fraction  of  the  cycle 
in  degrees  by  which  the  positive  peak  of  the  signal  leads  "zero  phase". 

This  "zero  phase"  is  arbitrarily  taken  as  occuring  exactly  at  half  the 
interval  between  two  consecutive  triggering  spikes.  At  this  instant  (it 
was  arranged  by  carefully  setting  the  five-tooth  gear  position,  relative  to 
the  propeller)  a propeller  blade  is  always  located  at  its  uppermost  posi- 
tion. The  way  phase  is  measured  is  better  illustrated  in  fig. (2). 

All  the  dynamic  test  reduced  data  is  given  in  graphical  form  in  the 
next  section  whilst  a listing  and  a sample  output  of  the  "dynamic  data  reduc 
tion  program"  is  given  in  Appendix  VI., 

Steady  State  Tests 

As  already  stated,  the  steady  state  tests  were  carried  out  more  as  a 
complement  to  the  dynamic  tests  than  to  investigate  in  detail  the  steady 
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state  characteristics  of  both  propeller  and  rudder. 

The  main  reason  behind  such  an  attitude  towards  steady  state  tests  was 
that  with  the  dynamometer  at  its  present  form,  there  was  not  much  more  that 
we  could  do.  A further  reduction  in  gain,  which  would  have  undoubtedly  in- 
creased our  testing  possibilities,  involved  tampering  with  the  electronics 
of  the  dynamometer  which  was  not  feasible  under  the  present  project. 

Thus,  as  soon  as  we  discovered  that  the  steady  loading  on  both  rudder 
and  propeller  was  unaffected  by  changes  in  clearance  the  steady  state  tests 
were  discontinued  and  we  concentrated  on  effects  on  dynamic  testing. 

The  testing  conditions  during  these  tests  were  also  quite  different  from 
the  ones  just  described.  First,  the  cavitation  number  being  so  low,  even  at 
the  smallest  amplifier  gain,  it  was  impossible  to  obtain  results  at  rudder 
angles  greater  than  10°.  Furthermore  the  amount  of  cavitation  appearing  on 
the  leading  edge  of  the  rudder  indicated  that  changes  in  the  tunnel  pressure 
would  have  no  measurable  effect  on  the  steady  state  loading.  Thus,  our  first 
decision  was  to  increase  the  pressure  in  the  tunnel  to  ’’atmospheric"  where  we 
knew  that  no  cavitation  would  occur  and  as  a result  the  dynamic  loads  on  the 
rudder  would  be  greatly  reduced. 

Sven  at  a pressure  as  high  as  this,  however,  we  were  still  unable  to  take 
any  readings  when  the  angle  of  attack  exceeded  20°,  without  saturating  the  ampli- 
fiers. As  a result,  at  such  angles,  in  order  to  decrease  the  rudder  loading, 
we  decided  to  reduce  the  values  of  both  n and  V^,  , keeping  their  ratios 
constant  (for  constant  J)  but  effectively  changing  once  again  the  value  of  O . 

The  second  testing  condition  that  was  altered  was  the  fact  that  no  wake 
screen  was  used.  There  are  two  main  reasons  behind  this  decision.  The  first 
one  is  that  at  the  time  we  carried  out  these  tests  there  was  no 
one  is  that  at  the  time  we  carried  out  these  tests  there  was  no 
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wake  screen  available;  the  prototype  had  failed  and  we  were  in  the  process 
of  constructing  the  second  (handwoven)  version  (see  Appendix  III).  The 
second  reason  was  that  we  felt  that  the  presence  of  a screen  would  not 
significantly  alter  the  influence  of  propeller  to  rudder  clearance  on  the 
steady  loading 

Rudder  Tests 

Had  it  not  been  for  the  significant  reduction  in  the  amplifier  gain, 
the  input  signals  used  in  the  steady  load  analysis  would  have  been  in  no 
way  different  from  the  ones  used  for  the  dynamic  tests.  In  fact,  if  we 
had  performed  a static  calibration  of  the  strain  gauges,  with  the  various 
amplifiers  set  at  their  maximum  gain,  nothing  (except  probably  lack  of 
time)  would  have  prevented  us  from  carrying  out  dynamic  and  steady  state 
tests  simultaneously. 

The  differences  in  the  two  tests  lie  in  the  way  these  signals  are 
analyzed  as  well  as  on  the  fact  that  for  steady  state  tests  we  examined 
the  output  of  five  rather  than  two  dynamometer  channels. 

The  instrument  used  to  measure  steady  loads  on  the  rudder  was  a di- 
gital indicator  (LEBOW  ASSOCIATES  model  7521).  This  instrument  was  origi- 
nally purchased  to  be  used  '.n  conjunction  with  the  load  cells  which  are  part 
of  another  steady  state  dynamometer  used  in  propeller  tunnel  as  well  as 
open  water  tests  at  M.I.T.  To  analyze  these  steady  state  loads  this  instru- 
ment incorporates  a transducer  unit  which  provides  the  excitation  volt- 
age to  the  load  cells  and  later  amplifies  and  filters  their  output  signals 
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prior  to  digitizing  them.  The  filter  and  the  digital  voltmeter  were  the 
only  parts  necessary  in  our  case  and  we  had  the  instrument  modified  by 
the  addition  of  another  unit  (essentially  a two-way  switch)  to  allow  us 
to  by-pass  the  transducer  circuit.  Thus,  the  indicator  could  be  used  as 
a digital  voltmeter  with  the  additional  feature  of  being  able  to  filter 
the  input  signal  prior  to  digitizing  it.  Thus,  by  allow* ng  a very  narrow 
band  of  frequencies  to  pass  through  the  filter,  the  fluctuation  of  the 
readings  would  be  kept  to  a minimum  and  we  had  no  problem  isolating 
the  steady  components  of  the  dynamometer  output  signals.  A high  degree 
of  filtering,  however,  required  a lot  of  damping  which  in  turn  meant  that 
a considerable  amount  of  time  had  to  be  allowed  before  being  able  to  take 
any  readings. 

There  was  one  additional  consideration  we  had  to  take  into  account 
before  using  this  instrument  for  the  analysis  of  the  signals.  This  is  the 
±1  volt  limitation  on  the  input  signal  voltage  imposed  by  the  instrument’s 
circuitry.  On  the  other  hand,  we  knew  that  the  output  from  the  amplifiers 
could  be  as  much  as  15  volts  (the  saturation  voltage  which  is  independent 
of  gain).  Thus  allowing  for  a safety  margin  we  introduced  a voltage  divider 
into  the  circuit  in  the  six-way  switch  capable  of  reducing  the  amplifier's 
output  voltage  to  1/17  of  the  original  value.  The  reduced  signal  was  then 
presented  tc  the  digital  indicator  at  the  "J^  " inlet  located  at  the  rear 
panel  of  the  instrument. 

Rudder  steady  state  tests  involved  measuring  the  output  of  five  out 
of  the  six  dynamometer  channels  for  thirteen  values  of  rudder  angle  (from 
+30°  to  -30°  at  5°  intervals)  and  at  four  longitudinal  clearances  between 
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propeller  and  rudder.  Only  one  lateral  position  was  examined  (the  design 
position).  The  output  of  the  sixth  dynamometer  channel  was  of  no  interest 
to  us  because  it  measured  vertical  forces  on  the  rudder  which  for  all  prac- 
tical purposes  are  negligible.  Prior  to  doing  any  testing  a static  cali- 
bration was  performed  to  evaluate  a matrix  of  influence  coefficients  re- 
lating the  applied  to  the  read  values  of  the  various  loads  (see  Appendix  IV). 

To  be  able  to  obtain  readings  at  large  angles  of  attack  the  no-load 
input  voltage  to  the  amplifiers  had  to  be  biased.  Because  of  the  reversal 
in  algebraic  sign  of  three  out  of  the  five  loads,  we  were  measuring, 

(Fz»Bx,  and  By),  as  we  moved  from  positive  to  negative  rudder  angles  it 
was  clear  that  the  same  no-load  bias  could  not  be  used  throughout  the  test 
without  saturating  the  amplifiers.  As  a result  the  rudder  steady  state 
tests  had  to  be  divided  into  two  groups,  one  dealing  with  positive  angles 
and  the  other  with  negative.  In  order  to  move  from  one  group  of  tests  to 
the  other  both  n and  had  to  be  turned  all  the  way  down  to  zero  and  the 
no-load  values  of  the  voltage  inputs  to  the  amplifiers  of  the  three  channels 
were  reversed. 

The  amount  of  bias  corresponding  to  each  amplifier  was  read  on  the 
digital  indicator  and  constituted  the  "zero  reading"  value  of  each  channel 
as  used  in  the  "steady  state  data  reduction  program"  (see  Appendix  VI). 

Propeller  Steady-State  Tests 

These  tests  were  carried  out  simultaneously  with  the  rudder  tests  and 
for  the  same  tunnel  conditions  and  number  of  points.  The  propeller  thrust 
and  torque  were  measured  by  two  different  load  cells  coming  as  one  package 


designed  and  built  to  be  used  in  this  propeller  tunnel.  Due  to  the  pro- 
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peller  rotation  the  load  cells  output  signal  is  transmitted  to  the  digital 
strain  gauge  indicators  (LEBOW  ASSOCIATES  model  66)  through  slip  rings. 
These  indicators  are  calibrated  so  that  they  measure  thrust  and  torque 
directly  so  that  no  further  calibration  was  needed.  A computer  program 
program  (not  included  here)  is  then  used  to  reduce  the  data  obtained 
giving  as  output  a set  of  ncn  dimensional  coefficients  including  Kp,  Kq, 
and  efficiency  n. 
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4.  EXPERIMENTAL  RESULTS 

Steady  State  Tests 

The  results  obtained  in  the  steady  rudder  experiments  were  very  much 
like  we  had  expected.  There  was  no  change  to  speak  of,  in  the  steady  state 
loading  at  any  of  the  four  longitudinal  clearances  we  investigated.  The 
lift  and  drag  coefficients  Cr  and  Cn  are  plotted  against  rudder  angle  at 
four  clearance  ratios.  The  four  curves  almost  completely  coincide  (see 
figs. (3,4  )).  As  already  stated  in  the  preceding  chapter  onl-  the  design 
lateral  location  of  rudder  was  investigated  (design  lateral  corresponds  to 
a displacement  ratio,  DR,  equal  to  .102). 

One  should  note  that  both  plots  indicate  quite  distinctly  that  zero 
lift,  and  hence  minimum  drag,  do  not  occur  at  exactly  0°  rudder  angle  but 
at  some  other  negative  angle.  This  is  not  unusual  but  actually  ouite  common 
to  all  ships  with  an  odd  number  of  propellers  or  any  number  of  unirotating 
propellers.  These  ships  in  order  to  be  able  to  sail  a straight  course 
must  set  their  rudder  at  some  angle  other  than  0°. 

Similar  conclusions  can  be  derived  from  the  propeller  steady  state 
tests.  Propeller  thrust  coefficient,  K , torque  coefficient,  K„,  and  ef- 
ficiency,n  , are  all  plotted  versus  clearance  ratio.  None  of  these 
quantities  appears  to  have  changed  to  any  observable  degree  with  propeller 
to  rudder  clearance,  (fig.(  5 )) 

We  can  thus  conclude  at  this  point  that  whatever  the  decision  on  the 
final  position  of  the  rudder,  it  will  have  to  be  made  on  the  basis  of  the 
results  of  the  dvnamic  tests  alone  since  no  constraints  appear  to  be  imposed 
by  steadv  state  loads. 
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Fig. 5:  Propeller  characteristics  versus  clearance  ratio.  Correct 
simulations  for  J for  ship  but  no  wake  screen  and  a is  too  high. 
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Dynamic  Test 

The  results  obtained  in  the  dynamic  tests  are  given  graphically  on 
figs.(  6-12).  The  first  four  plots,  one  at  each  lateral  location,  are  for 
the  same  tunnel  operating  conditions  and  in  all  of  them  the  model  to  ship 
dynamic  similarity  requirements  are  satisfied. 

The  air  content  of  the  water  is  kept  to  as  low  a value  as  possible  by  going 
through  the  procedure  described  previously  under  Dynamic  Test  Procedure.  It 
was,  however,  impossible  for  this  deaeration  procedure  to  have  exactly  the  same 
results  each  time,  so  that  reservations  to  compare  the  results  of  these  runs 
at  face  value  may  be  justifiable.  Nevertheless  we  did  try  to  keep  the  air- 
cor.tent  tc  a minimum  each  time  so  that  comparison  is  possible.  As  a result 
we  see  that  there  is  an  overall  drop  in  magnitude  of  the  loads  as  we  dis- 
place the  rudder  laterally  from  directly  behind  the  propeller  (DR=0). 


These  loads  reach  a minimum  somewhere  between  1"  and  2"  off  the  model  centre- 
line ( .102<DR<.201)  and  then  start  again  increasing  as  we  move  farther  to 


the  side. 


This  change  in  the  loading  may  be  interpreted  as  follows.  Exactly 


behind  the  propeller,  the  rudder  is  swept  by  the  propeller  huh  vortex 
and  the  induced  ve’.ocities  create  a fluctuating  pressure  field  which  in- 
creases the  vibratory  loading  on  the  rudder.  The  farther  away  we  are  from 
this  hub  vortex  the  smaller  its  effect  cr.  the  resulting  loads  until  we 
reach  a point  where  the  proximity  of  the  tip  vortex  envelope  causes  a simi- 
lar effect. 

When  we  consider  each  run  seoaratelv  we  see  that  the  orobabilitv  of 


the  air-content  of  the  water  chanping  within  its  duration  is  very  small. 
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Nevertheless,  we  decided  to  reduce  this  probability  even  sore  by  the  add- 
ition of  about  two  inches  of  water  to  the  top  of  the  tunnel,  thus  covering 
entirely  the  top  window.  This  was  done  because  earlier  tests  had  indicated 
that  this  window  was  the  major  source  of  air- leak  paths  into  the  tunnel - 

It  is  quite  clear  from  the  graphs  that  there  is  a substantial  variation 
of  both  the  force  and  the  moment  coefficients  (Kj.  and  Kg  respectively) 


with  changes  In  the  longitudinal  clearance  ratio,  €K- 
point  that  these  curves  agree  quite  well  with  those  of 
work.  [5,6] 


One  nay  note  at  this 
previous  experimental 


Close  inspection  of  the  curves  indicates  that  the  peaks  of  the  forces 
and  moments  occur  at  more  or  less  constant  intervals,  a fact  suggesting  that 
there  is  a definite  relationship  between  these  values  and  the  pitch  of  the 
tip  vortex  helix.  The  fact  that  all  four  curves  agree  in  that  respect 
strengthens  this  view.  The  pitch/no.  of  blades  of  the  helix  as  measured 
on  a photograph  (fig.  13)  is  approximately  5.1  cms,  which  corresponds  to  the 
size  of  the  interval  between  consecutive  peaks  or.  the  curves. 

In  other  words,  maximum  vibratory  loading  on  the  rudder  appears  to 
occur  whenever  the  tip  vortex  helix  "cuts"  the  rudder  at  a particular  point. 
Where  exactly  this  point  is  located  it  is  hard  to  tell.  Had  we  taken  a 
series  cf  photographs  of  the  tip  vortex  impinging  on  the  rudder  at  all 
the  clearance  ratios  It  might  have  probably  been  possible  to  determine  the 
exact  correlation  between  the  point  of  incidence  and  the  peaks  of  the  curves. 


Unfortunately,  this  was  not  done. 

There  is,  however,  an  additional  consideration  apparently  not  satisfied 


by  this  interpretation,  which  suggests  that  there  sight  be  an  additional 
source  of  excitation  of  these  loads,  this  is  the  fact  that  the  phase  of 
the  response  signal  appears  to  be  nearly  independent  of  the  distance,  between 
propeller  and  rudder  (see  fig. i-  ). 

Another  very  likely  source  of  vibration  Is  the  passage  of  a blade  be- 
hind the  wake  screen  five  tines  per  prop- Her  revolution  (i.e.  at  blade 
frequency).  This  scans  that  five  tiaes  per  revolution  the  loading  on  the 
propeller  changes  due  to  the  different  velocity  field  existing  behind  the 
screen.  As  a result  of  this,  a pressure  wave  is  generated  cen  each 
blade  which  travels  at  acaistic  speed  and  induces  vibratory  forces  oti  the 
rudder.  Because  of  such  high  speed,  the  few  cent  inters  between  the  various 
..c.nts  are  irssaterial  to  the  tine  it  takes  the  wave  to  reach  the  redder, 
fed  this  been  the  only  source  of  excitation  of  vibratory  the  rud- 

eer  then  one  would  not  expect  any  change  in  the  phase  of  the  response  signal. 
»it  we  have  already  interpreted  the  peaks  of  the  . »-ves  as  occuring 


ue 


the  tip  vortex  striking  the  rudder  a 


DoiflfS.  & 


quite  likely  that  the  total  loading  on  the  rudder  should  result  free  a oo»- 
w.nation  of  the  effects  of  the  two  sources.  The  net  load  on  the  rudder  wouJ 
then  be  a vector  sun  of  the  two  different  loads.  The  loads  generated  by 
no  presence  or  the  wake  screen,  being  predominant,  fixes  the  valise  of  the 
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occurs  when  the  tunnel  static  pressure  equals  the  barometric  pressure 
plus  the  hydrostatic  pressure  of  5 feet  of  water  (the  distance  from  the  top 
of  the  tunnel  to  the  centreline  of  the  test  section).  The  results  obtained 
at  this  pressure  are  shown  in  fig. ( 10).  From  this  graph  we  can  immediately 
observe  that  the  overall  shape  of  the  curve  has  changed  considerably,  there 
is  now  only  one  peak  much  less  pronounced  than  before,  and  that  the  mean 
level  of  vibration  has  dropped  to  about  1/3  of  its  value  at  low  pressure 
(c.f.  fig.  6). 

The  photographs  of  the  response  signal  also  show  another  difference. 

Now  the  phase  changes  with  clearance  [see  fig.  (11)].  It  actually  appears 
that  the  phase  goes  through  a complete  cycle  for  a change  in  longitudinal 
clearance  of  approximately  5cms  (equal  to  pitch/no.  of  blades  of  the  helix 
of  the  tip-vortex). 

One  can  ‘hus  conclude  that  at  increased  cavitation  numbers  the  vibratory 
loads  generated  by  the  blades  crossing  the  variable  wake  field  are  no 
longer  predominant  and  as  a result  of  that  the  contribution  of  the  tip-vortex 
induced  forces  has  increased. 

One  should  note  at  this  point,  that  even  though  it  is  not  visible 
at  high  pressure,  the  tip  vortex  is  always  present,  except  when  the  propeller 
is  delivering  no  thrust.  With  increasing  propeller  loading  the  vortex 
increases  in  intensity  and  its  core  starts  becoming  visible  when  the  pressure 
at  the  centre  is  reduced  until  it  reaches  vapour  pressure. 

An  extreme  case  of  the  effect  of  high  air-content  in  the  water  can 
be  observed  from  direct  comparison  of  fig.  (7)  and  (12),  Both  runs 


were  carried  out  under  the  same  tunnel  operating  conditions  but  in  the 
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later  case  the  air-content  was  so  high  that  the  propeller  was  barely  visible 
« in  the  tunnel.  Even  though  the  peaks  and  troughs  occur  at  approximately 

the  same  values  of  clearance  ratio  the  plots  show  considerable  differences 
in  the  overall  shape  of  the  response  curves. 

Further  insight  can  be  obtained  from  the  photographs  of  the  response 
signals.  There  we  noticed  that  the  repetitive  part  of  the  signal,  when 
the  air-content  was  high,  was  almost  completely  devoid  of  higher  order 
harmonics  of  blade  frequency  and  the  resulting  repetitive  waveform  was  al- 
most a pure  sinusoid.  This,  however,  was  not  the  case  when  the  air-content 
was  low,  suggesting  that  air  tends  to  dampen  higher  harmonics  and  that  the 
results  of  fig. (12)  are  essentially  the  pure  blade  frequency  components 
of  the  signal,  at  the  different  clearance  ratios.  We  do  not  know,  however, 
to  what  extent  this  is  true,  and  the  idea  of  pursuing  this  study  to  find 
out  what  happens  at  intermediate  values  of  air-content  was  almost  immediately 
abandoned  ■'  on  the  grounds  that  we  had  no  means  to  either  control  or  mea- 
sure the  quantity  of  air  in  the  tunnel.  As  a result  we  concentrated  our 
efforts  to  operate  at  as  low  an  air-content  as  possible  thus  obtaining  the 
response  of  the  dynamometer  to  all  the  components  of  the  propeller  excita- 
tion. 

Higher  order  harmonics  of  blade  frequency,  however,  introduce  an  ad- 
ditional problem.  Due  to  their  presence  the  response  signal  was  frequently 
distorted  to  an  extent  where  phase  was  not  always  easy  to  determine.  This 
was  particulary  the  case  at  clearance  ratios  corresponding  to  the  two  smallest 
clearance  ratios  (CR  equai  to  .33  and  .37)  and  the  two  largest  (CR  equal 
to  1.01  and  1.09).  In  all  these  cases  the  signal  was  considerably  distorted 
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to  the  point  that  the  location  of  the  positive  peak  (to  which  phase  was 
measured)  was  not  well  defined.  More  than  once  the  response  signal  had 
two  positive  peaks.  This  probably  explains  the  differences  in  the  shape 
of  the  four  curves  at  these  extreme  points . 

The  results  from  the  dynamic  tests  investigating  the  effect  of  changes 
in  rude  r angle  are  not  very  instructive.  Due  to  amplifier  saturation,  which 
was  caused  by  a considerable  increase  in  the  dynamic  loading  of  the  rudder, 
we  had  the  opportunity  of  examining  only  three  rudder  angles  (+5°  and  0°)  at 
two  clearance  ratios  (.57  and  .77).  At  both  positions  a minimum  occurred  at 

0°. 
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5 . CONCLUSIONS 

1.  The  results  of  these  experiments  indicate  clearly  that  propeller  induced 
vibratory  forces  are  strongly  influenced  by  propeller  cavitation.  The  growth 
and  collapse  of  cavitation  on  each  blade  as  its  passes  through  the  high  wake 
region  generates  a pressure  field  which  increases  the  rudder  vitratot'y  force 
by  a factor  of  three,  and  completely  alters  its  phase.  Realistic  design 
estimates  of  vibratory  force  levels  must  therefore  be  based  on  experiments 
or  theory  in  which  the  unsteady  cavitation  characteristics  of  the  propeller 
are  properly  accounted  for. 

2.  The  vibratory  force  on  the  rudder  is  critically  dependent  on  the  longi- 
tudinal clearance  between  the  propeller  and  rudder.  For  the  particular  design 
studied,  the  design  clearance  ratio  of  0.57  coincided  with  a local  maximum. 

A reduction  in  vibratory  force  by  a factor  of  ten  could  be  achieved  either  by 
reducing  the  clearance  ratio  to  0.41  or  increasing  it  to  0.68. 

3.  Transverse  clearance  between  the  propeller  and  rudder  also  influences 

the  vibratory  force  level.  In  this  case  the  design  value  of  DR =0.102  appears 
to  be  a good  choice.  Vibratory  forces  are  increased  both  for  smaller  and 
larger  values  of  the  transverse  clearance. 

4.  The  vibratory  force  on  the  rudder  at  zero  angle  of  attack  at  full  power 
has  an  amplitude  of  approximately  14,000  lbs.  at  the  design  location.  In 
the  absence  of  a structural  resonance , this  force  is  not  likely  to  cause 

r rob Lems  since  i.ls  magnitude  is  only  about  2%  of  the  maximum  steady  designed 
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force  on  the  rudder.  However,  any  structural  resonances  of  the  rudder  itself, 
of  the  steering  system  or  of  the  adjacent  hull  structure  at  blade  frequency 
or  twice  blade  frequency  could  result  in  vibration  problems. 

5.  Steady  state  propeller  and  rudder  characteristics  are  relatively  insensitive 
to  longitudinal  clearance  within  the  range  of  practical  possibility.  Conse- 
quently, a reduction  in  clearance  ratio  to  0.41  to  reduce  the  vibratory  force 
can  be  accomplished  without  reduction  in  propulsive  efficiency  or  steering 
effectiveness. 


6.  There  is  reason  to  believe  that  the  vibratory  force  on  the  rudder  may 
increase  by  an  order  of  magnitude  at  large  rudder  angles.  We  were  just  at 
the  point  of  acquiring  data  under  these  conditions  when  the  strain  gage  water- 
proofing failed.  Repairing  the  waterproofing  and  setting  the  experiment  up 

in  the  tunnel  was  unfortunately  not  feasible  within  the  existing  project. 

7.  The  shaft  inclination  of  the  actual  ship  was  not  modelled  in  this  experi- 
ment since  the  existing  propeller  drive  can  only  be  located  along  the  tunnel 
centerline.  It  was  anticipated  that  this  would  introduce  some  uncertainty 

in  interpreting  longitudinal  clearance.  However,  there  is  a further  problem 
associated  with  effect  of  shaft  inclination  on  unsteady  propeller  cavitation. 
Since  we  now  find  that  unsteady  propeller  cavitation  is  a major  variable, 
the  present  data,  run  with  zero  shaft  inclination,  may  underestimate  vibratory 


forces . 


fig.  a : Rudder  vibratory  load-.ng  versus  clearance  ratio,  DR=C.204 
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Rudder  vibratory  loading  versus  clearance  ratio,  DR 
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APPENDIX  I 

Testing  Facilities  and  Resulting  Design  Considerations 

The  various  tests  to  be  described  were  conducted  in  the  M.I.T.  variable 
pressure  water  tunnel.  The  tunnel  is  of  the  closed  circuit  type  having 
a test  cross-section  of  20  inches  square  and  is  capable  of  providing  water 
speeds  of  about  lOm/sec  (33ft/sec)  and  can  operate  at  pressures  varying 
from  atmospheric  (760mm  Hg)  to  about  100mm  Hg. 

All  this  information  is  relevant  since  it  provides  additional  restric- 
tions to  the  degree  of  simulation  which  can  be  achieved  in  the  experiment. 

Considerations  such  as  tunnel  wall  interference  and  longitudinal  travel 
of  the  propeller  drive  had  to  be  investigated  at  this  stage  and  the  fol- 
lowing decisions  were  made. 

The  transverse  positions  of  rudder  to  be  examined  were  to  be  restric- 
ted to  3 inches  each  way  of  the  centre  line  at  1 inch  intervals.  This 
would  provide  sufficient  points  to  enable  us  to  reach  acceptable  conclu- 
sions. This  decision  was  taken  to  satisfy  the  wall  interference  consider- 
ations as  well  as  practical  design  problems  which  would  have  arisen  other- 
wise. 

Longitudinally,  the  restriction  is  imposed  only  in  the  aft  movement 
of  the  rudder.  In  fact,  it  was  decided  that  it  would  be  much  more  simple 
to  move  the  propeller  longitudinally  relative  to  the  rudder,  rather  than 
the  other  way  around.  T5" ' s is  true  because  the  tunnel  is  built  to  make 
provision  for  longitudinal  movement  of  the  propeller.  Thus,  although  there 
is  a physical  limit  to  the  length  of  travel  of  the  drive  shaft  in  the  tunnel, 
it  is  such  that  it  did  not  interfere  with  any  of  the  positions  we  wanted 
to  examine.  Quarter  inch  spacing  near  the  design  location  and  half  inch 
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intervals  farther  away,  were  found  to  be  more  than  adequate. 

An  additional  provision  was  also  to  be  made  in  the  design  of  the 
rudder  mounting.  This  was  to  allow  for  variation  of  angle  of  attack  of 
the  rudder,  30°  each  way  at  5°  intervals. 

Finally,  we  were  to  try  to  test  at  as  high  a propeller  r.p.m.  as 
possible.  This  would  enable  us  to  operate  at  high  Reynold's  numbers  and, 
hence,  reduce  the  scale  effect  on  both  propeller  and  rudder.  It  has  also 
been  noticed  that  better  results  are  obtained  when  the  tunnel  is  operated 
at  higher  pressures. 

Natural  Frequency  Problems 

About  30  different  mounting  configurations  were  examined  in  order  to 
determine  which  one  provided  the  "best"  natural  frequency  characteristics 
on  the  rudder.  One  should  note,  however,  that  none  of  these  bore  all  the 
features  required  for  the  final  design,  since,  at  that  time,  we  were  only 
interested  in  dynamic  properties. 

The  experimental  determination  of  the  natural  frequency  of  each  con- 
figuration was  quite  simple.  Essentially,  the  method  used  was  the  same 
as  for  the  dynamic  calibration  (see  Appendix  IV)  with  the  difference  that 
the  frequency  of  excitation  was  varied  to  scan  the  entire  range  of  feasible 
operating  values  (0-300  Hz).  Since  we  were  not  interested  in  quantitative 
results,  the  output  signal  from  the  sensor  was  fed  into  an  ordinary  oscil- 
loscope where  amplitude  peaks  could  easily  be  seen. 

In  conjunction  with  the  experimental  investigation  of  the  system's 
dynamic  properties,  attempts  were  made  to  use  a theoretical  approach  to 
check  the  results.  Most  of  these  attempts,  however,  failed.  The  number 
of  approximations  as  well  as  the  fact  that  in  most  cases  the  boundary 
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conditions  were  ill  defined,  lead  to  results  which  were  seldom  verified 
experimentally . 

The  experience  gathered  from  the  various  tests  can  be  summarized  as 
follows : 

1)  The  natural  frequency  of  the  system  was  considerably  increased  when  the 
dynamometer  was  clamped  tight  at  more  than  one  location,  one  of  which  was 
as  close  to  the  rudder  as  possible. 

2)  In  most  cases  an  increase  in  the  mass  on  which  the  clamp  was  attached 
raised  the  natural  frequency  (not  considerably  however). 

3)  An  increase  in  damping  would  almost  always  result  in  a drop  in  the 
peak  response  of  the  system  although  this  would  also  increase  the  spread 
of  the  peak.  The  degree  to  which  the  natural  frequency  would  be  affected 
by  damping  was  very  much  dependent  on  the  overall  configuration. 

4)  The  window  boundary  conditions  can  have  little  or  great  effect  on  the 
dynamic  properties  of  the  entire  system  depending  on  the  way  the  dynamometer 
is  held. 

Unfortunately,  this  piece  of  valuable  information  was  not  perceived 
until  the  final  design  was  built.  Most  of  the  testing  until  then  had  been 
carried  out  with  the  window  sitting  on  a bench  rather  than  the  tunnel. 

This  was  decided  after  our  first  few  sets  of  results  obtained  on  the  tun- 
nel were  duplicated  on  the  bench. 

Bearing  in  mind  (1),  (2),  and  (3)  above  and  ignoring  (4)  we  went  on 
and  built  our  apparatus,  an  assembly  of  which  is  shown  in  fig.l.  The  effect 
of  the  window's  boundary  conditions  showed  up  Immediately  after  the  first 
test.  The  natural  frequency  of  the  window  just  resting  on  top  of  the  tunnel 
was  210  Hz  (not  far  from  what  we  were  e):pectlng). 
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However,  as  soon  as  the  window  was  bolted  down  tight,  the  same 
natural  frequency  dropped  to  93Hz.  Unfortunately,  to  avoid  leakage  through 
the  window  sides,  it  was  in  this  latter  way  that  the  experiment  had  to  be 
carried  out.  One  should  note  that  both  values  quoted  above  were  for  tests 
in  air.  In  water  the"added  mass  effect"  reduced  these  values  even  further. 

This  sudden  drop  in  the  first  natural  frequency,  it  was  later  decided, 
was  due  to  the  entire  rig  vibrating  on  top  of  the  window  about  a horizontal 
axis  passing  through  the  centre  of  the  plate  (see  fig.l).  This  mode  of 
vibration  was  present  only  when  the  window  was  constrained  so  as  not  to 
move  on  the  tunnel.  The  second  natural  frequency,  however,  was  due  to  the 
2-noded  flexural  vibration  of  the  sensor-rudder  subassembly  acting  as  a 
very  flexible  cantilevered  beam  carrying  a heavy  mass  at  its  end. 

Thus  our  original  aim  of  operating  at  a frequency  below  the  first 
natural  had  to  be  abandoned  and  we  had  to  live  with  the  fact  that  we  had 
to  operate  between  the  first  and  second  modes. 

Operating  between  two  natural  frequencies  is  frequently  avoided  be- 
cause of  instabilities  in  amplitude  and  mainly  phase  which  may  occur  in 
the  output  response.  Having  established,  however,  that  the  first  natural 
frequency  was  not  due  to  the  rudder-sensor  subassembly  but  rather  due  to 
the  entire  rig,  we  anticipated  little  variation  in  either  amplitude  or 
phase  provided  the  distance  between  the  two  peaks  was  sufficiently  large. 

To  ensure  an  even  larger  spacing,  we  decided  to  lower  the  stiffness 
of  the  system  by  soft  mounting  the  steel  plates  onto  the  window. 

This  decision  was  proved  correct.  When  we  scanned  the  range  of  pos- 
sible operating  frequencies  we  observed  that  the  first  natural  frequency 
had  dropped  to  76Hz  (in  water)  whilst  the  second  had  been  pushed  up  to  275  Hz, 


- cn  _ 


We  also  discovered  that  between  120  and  200  Hz  there  existed  a band  in 
which  neither  amplitude  nor  phase  appeared  to  change  at  all. 

Additional  tests,  using  a ubiquitous  spectrum  analyzer  (FEDERAL  SCIEN- 
TIFIC MODEL  VA-15A),  to  scan  the  frequencies  between  these  two  values  made 
us  decide  to  fix  our  operating  blade  frequency  for  the  experiment  at  125  Hz 
which  corresponds  to  a shaft  frequency,  N,  equal  to  1500  R.P.K. 

Having  fixed  one  of  the  tunnel  operating  parameters,  dynamic  simila- 
rity considerations  immediately  determine  the  remaining  two,  namely  water 

soeed,  V , and  static  pressure,  P (see  AnDendix  V ).  From  fig. (22) 

T ST 

and  ( 23  ) we  see  that  the  resulting  final  values  of  these  two  parameters 
are  17.3  ft/sec  and  236mmHg  respectively. 
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APPENDIX  II 

The  Six  Component  Dynamometer 

The  six  degrees  of  freedom  dynamometer  was  designed  and  built  between 
1369  and  1970,  as  part  of  a project  sponsored  by  the  National  Science 
Foundation.  Its  original  purpose  was  to  measure  unsteady  forces  and  moments 
on  c.  propeller.  A description  of  the  sensor  and  the  strain  guages'  mounting 
and  wiring  is  given  in  great  detail  in  refs.  (11-13)  and  is  beyond  the  scone 
of  this  report. 

It  is,  however,  necessary  to  outiire  the  few  alterations  on  the  original 
design  that  had  to  be  made  to  enable  its  use  as  a rudder  dynamometer  capable 
of  measuring  both  steady  and  vibratory  loads.  A couple  more  points  relevant 
to  this  experiment  are  also  brought  up  for  the  benefit  of  similar  tests  in 
the  future. 

The  most  important  modification  on  the  original  design  was  the  intro- 
duction of  potentiometers  in  the  amplifier  circuits.  These  are  essentially, 
additional  sources  of  input  to  the  amplifiers  which  can  be  controlled  ex- 
ternally and  can  be  used  to  alter  the  net  amplifier  output  voltage.  With 
their  addition,  the  no-load  input  voltage  can  be  biased  in  such  a way  so 
that  when  loads  are  applied,  the  net  output  voltage  will  be  below  the  ampli- 
fier's saturation  value.  Therefore,  depending  on  the  value  of  the  potentio- 
meter resistors,  the  span  of  the  possible  voltages  which  may  be  obtained 
from  the  amplifier,  can  be  increased  to  as  much  as  twice  the  value  without 
the  potentiometers. 


At  present,  the  potentiometers  for  each  strain  gauge  bridge  are 
capable  of  causing,  at  minimum  gain,  the  bias  shown  on  the  table  below. 
Clearly,  voltages  greater  than  15  + |bias|  will  always  saturate  the  amp- 
lifier (this  was  frequently  the  case  when  the  rudder  angle  of  attack  was 
greater  than  5°);  it  should  also  be  noted  that  the  higher  the  amplifier 
gain  the  smaller  the  voltage  needed  to  cause  saturation.  A schematic  of 
the  strain  gauge  bridge  circuit  is  shown  in  fig.(  15). 
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The  amount  of  bias  which  can  be  achieved  by  the  individual  potentio- 
meters can  be  further  increased  by  reducing  the  value  of  R^.  There  is, 
however,  no  easy  access  to  this  resistor  so  that  whenever  possible  such  a 
change  should  Le  left  as  a last  resort. 

Lowering  the  amplifier  gain  is  another  way  of  avoiding  saturation. 

This  is  a function  of  the  ratio  R^g/Rge  an<3  the  feedback  network  resistor 
Rr.  A reduction  in  either  of  these  quantities  will  result  in  lowering 
the  gain.  The  easier  one  to  reduce  is  RAg  by  simply  adding  another  resistor 
R across  the  terminals  1 and  2 (see  fig.  15  ).  The  minimum  gain  which  can 
be  achieved  in  this  way  is  1/10  of  the  original  value.  Further  reductions 
in  gain  are  much  harder  to  accomplish  and  require,  a major  overhaul  of  the 
dynamometer. 

To  prevent  damage  to  the  circuitry  the  strain  gauges  had  to  be  care- 
fully waterproofed  (previous  experiments,  using  this  dynamometer,  had  failed 


dynamometer  strain  gauge  circuit 


because  of  unsuccessful  waterproofing).  This  was  accomplished  by  lowering 
the  bottom  part  of  the  sensor  into  molten  bee’s  wax  until  a thick  protective 
layer  was  built  up  on  the  gauges.  Excess  wax  was  machined  off  on  a lathe. 
The  wax  itself  was  coated  with  a layer  of  M-COAT  G resin  (M-LIHE  ACCESSORIES) 
which  although  it  proved  out  to  adhere  poorly  on  either  wax  or  steel,  it 
nevertheless  provided  some  additional  shock  protection  of  the  gauges. 

The  dynamometer  before  and  after  waterproofing  can  be  seen  In  f Ig.  ( 16-  ) - 

Some  additional  wax  was  used  to  fill  part  of  the  core  of  the  dynamo- 
meter. This  was  done  in  appreciation  of  the  fact  that  the  reduced  pressure 
in  the  tunnel  could  force  air  bubbles  through  the  small  holes  on  the  dyna- 
mometer wall  (in  the  proximity  of  the  gauges)  into  the  protective  layer 
of  wax,  thus  damaging  the  waterproofing. 

Finally,  since  the  dynamometer  was  now  used  to  measure  forces  and 
moments  on  a rudder  the  propeller  based  coordinate  system  was  no  longer 
valid.  As  a result,  the  various  strain  gauge  bridges,  channels,  etc.  had 
to  be  relabelled  in  a manner  consistent  with  previous  rudder  tests.  The 


"new"  coordinate  system  adopted  is  the  universally  accepted 


handed  system  (see  fig.(  17)). 


ie  Waveform  Eductor 


time  the  repetitive  waveform  occurs,  separately  integrating  to  suppress 
the  noise,  and  finally  storing  each  integrated  sample  in  its  memory.  With 
each  sweep  the  contents  of  the  memory  are  reviewed  and  asymptotically 
approach  the  point  by  point  value  of  the  repetitive  part  of  the  signal  even 
though  at  the  input  the  repetitive  part  of  the  signal  might  have  been  com- 
pletely buried  in  noise. 

Synchronizing  the  sweep  with  the  external  event  can  be  achieved  by 
supplying  a trigger  pulse  to  the  Eductor  either  automatically  (externally 
or  internally)  or  manually.  (Since  in  the  experiment  to  be  described  ex- 
ternal triggering  was  used  there  is  no  need  at  this  point  to  elaborate  any 
further  on  the  other  triggering  methods.) 

An  important  feature  of  the  Eductor  is  that  it  allows  the  user  to  . 
choose  both  the  sweep  duration  and  the  instrument's  characteristic  time 
constant.  In  fact  these  quantities  are  interrelated  and  it  is  possible 
by  combining  them  in  a proper  way  to  optimize  the  efficiency  of  a run. 

The  importance  of  this  feature  is  evident  when  one  considers  the  number 
of  points  that  had  to  be  investigated  in  each  run. 

By  sweep  duration  one  means  that  fraction  of  the  sweep  period  (or 
time  between  two  consecutive  triggering  spikes)  the  machine  examines  the 
input  signal  and  adds  information  to  its  memory.  For  efficient  runs  one 
should  arrange  so  that  the  sweep  duration  is  just  less  than  the  sweep  per- 
iod. The  reason  behind  this  can  best  be  illustrated  by  considering  the 
following 

For  a 5-bladed  propeller  the  period  of  the  repetitive  waveform  is 
given  by  (1/blade  frequency)  and  is  essentially  the  time  it  takes  for  one 
blade  to  cover  1/5  of  a revolution  (i.e.  72°).  For  a blade  frequency  of 
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125Hz  this  period  is  thus  8msecs.  Assuming  r.ow  that  we  are  interested  in 
viewing  only  one  period  of  this  wave,  (i.e.  we  want  the  sweep  period  to 
be  equal  to  the  waveform  period).  If  the  swee^  duration  is  set  to  7.8msecs, 
say,  the  Eductor  will  be  adding  to  its  memory  over  97%  of  the  signal  dura- 
tion having  to  wait  for  only  2msecs  for  the  next  trigger.  On  the  scope 
the  signal  will  look  as  on  fig.  (18)  case  a_. 

Had  the  sweep  duration  been  a fraction  larger  than  the  waveform  period, 
8.2msecs  say,  then  the  Eductor  would  miss  completely  one  triggering  spike 
taking,  thus,  twice  as  much  time  to  extract  the  complete  signal.  On  the 
scope  the  signal  would  now  look  as  on  fig.  (18)  case  b_.  Another  reason  why 
the  sweep  duration  in  case  a_  should  be  preferred  to  the  one  in  case  ]b  is 
strictly  qualitative.  It  is  essentially  due  to  the  fact  that  during  "off" 
time  there  is  a tendency  for  significant  quantities  of  charge  to  leak  off 
the  memory  capacitors  and  correspondingly  a greater  likelihood  of  distor- 
tion in  the  stored  waveform.  It  is  evident  that  in  the  second  case  the 
"off"  time  of  the  instrument  is  greater. 

The  ratio  (sweep  duration/sweep  period)  is  known  as  the  duty  factor 
(D.F.)  and  its  knowledge  is  important  in  estimating  the  real  time  required 
to  charge  the  memory  to  within  1%  of  the  applied  voltage. 

The  characteristic  time  constant  (t)  determines  directly  the  time 
required  for  the  Eductor  to  store  the  repetitive  waveform  component  of 
the  input  signal  in  its  memory.  One  should  note,  however,  that  this  time 
is  actual  instrument  integrating  (or  sweeping)  time  and  differs  from  the 
!<-,il  or  "clock  i tine"  by  the  duty  f.iclop. 

The  i n*;t  rumen  t ' *:•  active  time  required  to  charge  1 he  memory  If'  1%  nl 
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the  applied  voltage  is  5*t.  Therefore,  real  time  is  given  by 


5*t 

Real  Time  = „ „ 


In  choosing  the  operating  value  of  t one  should  give  equal  considera- 
tion to  both  the  efficiency  of  the  run  and  to  the  fact  that  increasing 
t gives  a slightly  better  resolution.  For  example , for  a signal  with  large 
noise  to  signal  ratio  a larger  value  of  t should  be  preferred  even  though 
this  may  penalize  the  efficiency.  An  idea  of  how  to  pick  an  "optimum" 
value  of  t may  be  taken  from  the  instrument's  operating  manual  where  t 
is  plotted  vs.  D.F.,  sweep  duration  and  noise  to  signal  ratio. 

Measurement  of  the  Model  Testing  Parameters 

Having  already  decided  on  the  final  values  of  the  operating  parameters 

N,  V , and  P it  is  worth  mentioning  the  way  they  were  measured  in  the 
T st 

tunnel. 

Shaft  frequency  is  measured  by  means  of  a magnetic  pick-up  sending 
out  pulses,  everytime  a tooth  of  a 60-tooth  gear  passes  underneath  it. 

This  gear  which  rotates  with  a shaft,  comes  as  part  of  the  propeller  thrust 
end  torque  transducer  package,  located  at  the  rear  end  of  the  tunnel.  The 
number  of  pulses  thus  generated  is  measured  by  an  electronic  counter  (HP 
model  5321B)  over  a given  period  of  time  which  for  the  particular  case 
was  set  to  one  second.  Clearly  th>’  number  of  pulses  per  second  received 
by  the  counter  corresponds  to  the  shaft  frequency  measured  in  r.p.m. 

Water  speed  can  be  determined  from  the  difference  in  the  pressures 
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obtained  at  two  points  located  at  the  entry  to  the  test  section.  This 
pressure  difference,  in  this  particular  test,  was  measured  on  a 3 vard  long 
differential  manometer  which  uses  "blue  fluid"  as  the  working  fluid  instead 
of  mercury.  This  last  decision  was  made  in  anticipation  of  the  fact  that, 
in  order  to  increase  the  sensitivity  of  the  manometer  readings,  a low  spe- 
cific gravity  fluid  shoulG  be  used.  The  fluid  chosen  (MERIAM  D-8325  indi- 
cator fluid)  is  considerably  less  dense  than  mercury,  having  a specific 
gravity  of  1.75  at  55°F.  The  relationship  between  the  height  of  the  "blue 

fluid"  column,  H , and  the  water  soeed  at  the  entry  to  the  test  section, 
b 

V,  is  given  by 


H = 4.11  x cf  x V2-026  (1) 

b 

where  C'  is  a temperature  correction  factor. 

In  our  case,  however,  the  desired  water  speed  jf  17.3  ft/sec  is  the 
value  at  the  middle  of  the  test  section  which  due  to  the  presence  of  the 
splitter  plate  and  the  resulting  decrease  in  cross  sectional  area,  is  dif- 
ferent from  the  value  at  the  entry.  Measurements  with  a pitot  static  probe 
indicated  that  the  former  value  was  1.091  x V so  that  equation  (1)  above 
can  be  rewritten  as 


= u.H  x c'  x (V^/1.091)2'026 


(2) 


Substituting  17.31  ft/sec  for  V and  1.35264  for  C*  we  get  from  (2) 


that  the  value  of  H corresDonds  to  the  soeed  in  the  middle  of  the  test 
b 
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section  is  1506mm. 

Finally,  static  pressure  at  the  centreline  of  the  propeller  is  measured 
at  two  points,  one  at  the  entry  and  the  other  one  at  the  exit  of  the  test 
section.  The  two  pressures  transmitted  through  tubes  are  averaged  at  a 
merging  point,  and  the  resulting  mean  value  is  read  on  a mercury  manometer 
(reading  absolute  rather  than  gauge  pressure). 


APPENDIX  III 


Water  Tunnel  Wake  Survey 

To  achieve  an  even  higher  degree  of  simulation  it  was  decided  to 
reproduce  in  the  tunnel  the  variable  wake  obtained  behind  the  ship.  Such 
a wake  can  be  generated  by  inserting  at  the  entrance  to  the  test  section 
a screen  of  variable  solidity  capable  of  creating  a velocity  profile  having 
a variation  in  the  axial  velocity  component. 

The  wake  behind  the  ship  was  investigated  in  a wake  survey  carried  out 
on  a model  in  a towing  tank.  The  results  of  this  wake  survey  were  given 
to  us  as  ratios  of  the  axial  velocity  of  the  flow  behind  the  ship  to  the 
average  free  stream  velocity,  plotted  against  angular  position  at  4 different 
radial  distances  from  the  propeller  shaft  centreline. 

From  these  results  we  were  able  to  make  the  following  observations: 

1)  A significant  reduction  of  the  axial  velocity  component  was  evident  at 
all  radii  and  at  angles  varying  from  +20  to  -20  degrees  to  the  upward  ver- 
tical. This  reduction,  caused  by  potential  flow  pressure  gradient  inter- 
acting with  hull  viscosity  effects  as  well  as  the  pressure  of  the  propeller 
shaft  the  axis  of  which  was  inclined  relative  to  the  flow,  was  more  pro- 
nounced at  the  smaller  radii. 

2)  At  all  radii  there  existed  a transition  zone  between  angles  of  ±20-35° 
where  the  flow  was  evidently  affected  by  the  presence  of  the  vee-strut. 

3)  At  the  largest  radius,  1. 174xpropeller  radius  the  flow  did  not  appear 
to  recover  from  this  transition  zone,  until  the  angle  reached  ±60°.  This 
fact  and  close  inspection  of  the  ship's  drawings  suggested  that  at  such  a 


large  radius  the  proximity  of  the  hull  is  a critical  factor. 

4)  At  large  angles  the  longitudinal  velocity-component  ratio  appeared  to 
be  fairly  constant  for  all  radii. 

In  the  tunnel  matching  the  wake  of  the  ship  presented  certain  problems. 
Bared  upon  the  observations  above,  it  was  decided  to  build  a screen  which 
would  essentially  consist  of  a frame  of  trapezoidal  shape  (see  fig.20  ) 
on  which  a wire  mesh  whose  variable  solidity  was  to  be  evaluated  experi- 
mentally was  welded  (solidity  being  defined  as  blocked  area,  Ab,  divided 
by  total  area,  Ab  + Ao,  where  Ab  is  the  projected  area  of  the  screen  elements 
and  Ao  that  of  the  openings). 

We  thus  started  with  a screen  of  given  solidity  (.272)  and  carried 
out  a wake  survey  in  the  tunnel  in  which  we  measured  the  velocity  distribu- 
tion behind  the  screen  at  four  radii  (corresponding  to  the  four  radii  used  in 
the  towing  tank  survey).  The  water  speed  was  measured  by  means  of  a pitot  static 
probe,  attached  to  the  propeller  drive  shaft  and  thus  able  to  assume  any 
angular  position  by  merely  moving  the  drive  belts  at  the  back  of  the  tunnel 
(see  fig.19  ).  The  angle  obtained  was  read  to  the  nearest  degree  on  a pro- 
tractor fixed  at  the  end  of  the  shaft. 

At  each  radius  the  water  speed  was  measured  at  several  angles  within 
the  region  influenced  by  the  screen  and  the  ratio  of  each  of  these  values 
to  the  fully  developed  "free"  stream  velocity  was  computed.  Finally,  the 
average  of  all  these  ratios,  at  each  radius,  was  obtained  and  compared  to 
the  one  obtained  in  the  towing  tank  survey. 

Depending  on  the  difference  between  the  2 averages,  the  solidity  of 
the  screen  was  either  decreased  by  removing  some  of  the  wire,  or  increased 
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by  adding  bits  of  wire  mesh  wherever  necessary.  This  procedure  was  repeated 
until  obtained  and  desired  velocity  distributions  bore  satisfactorily  close 
resemblance . 

The  drawbacks  of  this  method  became  apparent  during  the  first  test 
runs.  Satisfying  the  high  velocity  requirements  at  certain  radii  meant 
reducing  the  solidity  of  the  screen  to  a point  where  many  and  sometimes 
all  of  the  transverse  members  of  the  mesh  had  to  be  removed.  Besides  weak- 
ening the  screen  as  a whole,  this  meant  that  the  remaining  longitudinal 
wires  were  now  free  to  deflect  and  reach  larger  curvatures  than  before 
with  the  result  of  being  subjected  to  higher  stresses. 

This  combined  with  the  increase  in  the  loading  due  to  the  reduction 
in  pressure  was  sufficient  to  cause  fatigue  failure  in  a number  of  screen 
wires. 

All  this  prompted  us  to  find  some  other  way  of  obtaining  the  required 
variable  solidity  and  so  we  finally  decided  to  handweave  a screen  on  the 
frame . To  do  so  we  drilled  a large  number  of  closely  spaced  holes  on  all 
four  sides  of  the  frame  through  which  long  pieces  of  wire  were  weaved 
forming  a mesh  of  longitudinal  and  transverse  wires  which  now  had  consider- 
ably greater  strength. 

The  first  version  of  the  screen  was  deliberately  made  less  solid  than 
what  we  anticipated  the  final  version  (after  the  wake  survey)  would  be,  to 
make  sure  that  any  changes  would  involve  additions  rather  than  removal  of 
wire.  The  final  version  of  the  screen  is  shown  on  fig.  (20). 

The  resulting  wake  behind  the  final  screen  bore  a satisfactory  resem- 
blance with  the  wake  obtained  behind  the  model  in  the  towing  tank  survey. 

The  results  of  the  two  surveys  are  given  in  the  following  t^ble  (I). 


. n wake  screen  and  (b)  complete  model 
Fig. 20:  (a)  Handwoven  wake  scree 

scale  simulation  of  ship  configuration. 
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Splitter  Plate  Boundary  Layer 


The  introduction  of  the  splitter  plate,  as  already  stated,  was  decided 
with  2 considerations  in  mind.  Firstly,  it  was  to  serve  the  purpose  of 
simulating  the  ship’s  hull  by  providing  the  proper  propeller  tip  clear- 
ance. Secondly,  it  was  hoped  that  by  making  it  extend  less  upstream  than 
the  test  section  wall,  the  resulting  boundary  layer  would  be  of  smaller 
thickness  than  the  one  that  would  have  developed  on  the  wall  otherwise. 

For  this,  however,  to  be  the  case,  it  was  established  in  previous 
experiments,  that  the  leading  edge  had  to  be  slightly  curved  upwards.  The 
decree  of  this  inclination  was  decided  empirically  based  on  the  experience 
obtained  with  an  earlier  splitter  plate  used  in  the  tunnel  f*)j,  whose 
30°  leading  edge  angle  had  managed  to  create  a boundary  layer  of  relatively 
small  thickness. 

Measuring  the  boundary  layer  thickness  on  the  tunnel  was  relat ively 
simple,  involving  a method  quite  similar  to  that  used  in  the  wuxe  survey. 
The  pitot  static  probe  was  raised  until  it  almost  touched  the  splitter 
plate  at  0°  to  the  vertical..  The  water  speed  was  measured  at  this  radius 
and,  as  it  was  expected,  was  found  to  be  very  small  . The  probe  was  then 
rotated  by  small  angular  movements  until  a point  was  reached  from  v.T  ;'  r. 
or.  the  velocity  read  was  constant.  This  point  marked  the  end  of  the  bound - 
trv  "r.^r.  ~'vy  polar  coordinates  (r,€-)  of  this  rcint  one  ran  evaluate 
the  boundary  layer  thickness  using  the  expression 

6 = r(l-cos9)  + h 


where  h is  the  1/1 & of  an  inch  clearance  between  the  probe  at  its  upper- 
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cost  position  and  the  splitter  p'ate. 

For  r = 7.5625  inches  and  9 = 16°  the  boundary  layer  thickness  was 
found  to  be  .355  inches  which  is  indeed  a very  satisfactory  value.  This 
result  also  happens  to  agree  remarkably  well  with  the  one  obtained  if  the 
following  expression  for  a turbulent  boundary  layer  on  a flat  plate  is 
used 

. _ 0.37x 

6 ‘ 1/5 

R 


where  p-  Vd  is  the  Reynolds ' nusber 
v 


x=d  is  the  distance  of  the  axis  of  rotation  of  the  rudder  frees 
the  leading  edge  of  the  splitter  plate.  (Ft) 

V is  water  speed  (ft/sec). 

v is  ksneratic  viscosity  of  water  at  21°C  ft  /sec. 


Substituting  then  in  the  above  equation  V = 19.76  ft/sec,  x=d=!9/12  ft 
and  v = 1.059xio  ^ , we  get  that  the  boundary  layer  thickness  on  the  plate 


should  be  .357  inches. 
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APPENDIX  IV 


Dynamic  Calibration 

As  already  stated,  a dynamic  calibration  was  carried  out  to  determine 
the  complex  coefficients  relating  the  true  vibratory  force  and  the  moment 
(F  and  Bv  respectively)  on  the  rudder  to  the  output  of  the  corresponding 
strain  gauges. 

A desired  calibration  force  was  produced  by  a known  alternating  cur- 
rent in  a coil  placed  in  the  air  gap  of  a permanent  horseshoe  shaped  magnet 
and  attached  to  the  rudder  at  given  locations. 

The  coil-magnet  combination  was  first  calibrated  by  measuring  the 
steady  force  produced  by  a known  direct  current  in  the  coil.  The  constant 
of  the  particular  coil  we  used  had  been  found  by  F.M.  Lewis  [6j  to  be 
.775  lbs/amp.  Recalibration  using  a chemical  scale  with  weights  balancing 
out  the  steady  magnetic  force,  yielded  a constant  equal  to  .75  lbs/amp. 

The  difference  between  the  two  values  was  attributed  to  the  fact  that  the 
magnet  had  been  left  for  lengthy  periods  without  its  keeper  with  the  result 
of  losing  part  of  the  strength  of  its  magnetic  field. 

The  strain  gauge  calibration  procedure  was  as  follows.  The  coil  was 
first  bolted  onto  the  rudder  at  some  point  Y^  inches  from  the  splitter 
plate  (which  was  used  as  datum)  as  shown  in  fig.  ( 21).  With  the  magnet 
in  position  an  alternating  current  from  a signal  generator  (HP  model  204C) 
of  known  frequency  was  passed  through  the  coil  via  an  amplifier  (MC  INTOSH 
model  200).  The  amplitude  of  the  field  current  was  determined  by  r ensuring 
V,  the  R.M.S.  voltage  drop  across  a 4.6ft  resistor  placed  in  series  with  the 


coil.  (V  was  measured  with  an  HP  vacuum  tube  voltmeter  model  400D).  Such 
an  arrangement  allowed  the  operator  to  keep  the  current  input  to  the  coil 
at  a constant  value  independent  of  the  exciting  frequency. 

The  resulting  peak  force  on  the  rudder  was  thus: 


r i /;  * .75  xv 

CAL o lbs- 

The  response  of  the  strain  gauges,  F and  to  this  force  was 

then  fed  through  the  waveform  eductor  onto  an  oscilloscope.  A Polaroid 
picture  of  this  output  signal  was  then  taken,  on  which  the  amplitude  and 
phase  were  read. 

This  procedure  was  then  repealed  for  another  position  of  the  coil,  Y^ 
inches  from  the  splitter  plate.  Given  the  amplitude  and  phase  of  the  sensor 
responses  the  complex  coefficients , ai j , can  be  evaluated  from  the  following 
equations ; 
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The  solution  of  these  four  equations,  which  involves  complex  algebra, 
is  carried  out  in  the  first  part  of  the  "dynamic  data  reduction  program" 
(see  Appendix  VI). 

Static  Calibration 


A static  calibration  of  the  sensor  gauges  was  carried  out  to  determine 
their  response  to  steady  loads.  Although  in  theory  it  could  have  been  carried 


74 


out  in  much  the  same  way  as  the  dynamic  calibration  (by  passing  a direct 
current  through  the  coil)  a different  method  was  used.  This  was  done  in 
anticipation  of  the  fact  that  during  actual  testing  the  steady  loads  acting 
on  the  rudder  would  be  of  a much  greater  order  of  magnitude  than  the  ones 
that  would  possibly  be  applied  through  the  coil. 

The  method  used  essentially  consisted  of  hanging  weights  from  a speci- 
ally designed  calibration  rig  bolted  to  the  dynamometer  in  place  of  the 
rudder  (fig. (21)).  The  rig  makes  provision  for  hanging  weights  from  a 
number  of  different  points  at  various  distances  from  the  gauges  so  that, 
depending  on  the  point  of  application,  a net  force  or  a net  moment  may  be 
applied.  Since  it  was  already  decided  that  the  vertically  upward  force 
on  the  rudder,  F , was  not  to  be  measured,  the  corresponding  channel  was 
not  calibrated. 

Having  a fairly  good  idea  of  the  expected  maximum  values  of  the  various 
loads  we  adjusted  the  amplifiers'  gain  accordingly  by  introducing  resis- 
tances, Rg,  across  the  terminals  1 and  2 of  each  bridge  circuit  (see  fig. 15). 
The  values  of  these  resistors  and  the  resulting  reduction  in  gain  (as  a 
fraction  of  maximum  possible  gain)  is  given  below. 
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With  the  pain  of  each  remaining  channel  adjusted  to  these  final  values 
the  response  of  the  strain  gauges  to  various  loads  was  read  on  a digital 
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indicator  (LEBOW  ASSOC,  model  7521).  As  in  the  case  of  the  dynamic  calib- 
ration the  net  force  (or  moment)  on  the  dynamometer  in  one  direction  was 
an  algebraic  combination  of  the  readings  of  the  force  and  moment  channels 
in  the  same  direction. 

The  results  of  the  calibration,  namely  the  matrix  of  the  coefficients 
relating  actual  to  read  values,  are  shown  below. 


F 

X 

-0.077 

0 

0 

-0.063 

0 

F | 

XT 

F 

Z 

0 

0.072 

0.137 

0 

0 

! 

F t 

zr  j 

B 

X 

— 

0 

-0.051 

0.455 

0 

0 

B j 

XT  j 

B 

z 

-0.043 

0 

0 

0.223 

0 

B 

zr 

B 

_ y _ 

0 

0 

0 

0 

-.65 

B 

L ^ J 

These  results  are  used  in  the  "steady  state  data  reduction  program" 
in  the  evaluation  of  the  different  forces  and  moments  acting  on  the 
rudder  (see  Appendix  VI). 
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APPENDIX  V 


Model  to  Ship  Dynamic  Similarity  Considerations 

In  order  for  the  results  of  the  model  experiments  to  be  of  any  value 
to  the  actual  ship,  there  are  certain  laws  governing  model  to  ship  simili- 
tude which  need  to  be  obeyed.  These  laws  are  determined  by  making  use  of 
dimensional  analysis. 

If  we  apply  dimensional  analysis  in  examining  the  thrust  on  a propeller 
we  come  up  with  an  equation  for  the  thrust  coefficient  of  the  form 


2 

PD  V. 


=f  r_j5. , 

> [v  v 


nD 

V. 


P 

pv. 


21  = f[Fn.  V J.  » ] 


This  equation  states  that  for  similar  flow  patterns  and  thrust  coeffi- 
cients for  the  model  and  the  ship,  the  values  of  the  non-dimensional  groups 
on  the  right-hand  side  in  the  two  cases  must  be  the  same. 

Satisfying  both  the  Froude  number  and  the  Reynolds  number  equality 
constraints  is  not  possible.  In  water  tunnel  experiments,  (as  opposed  to 
open  water  tests)  it  is  common  practice  to  satisfy  the  Reynolds  number  re- 
quirement, being  prepared  to  accept  differences  in  the  local  pressure  dis- 
tributions across  the  disk,  for  the  benefit  of  avoiding  more  critical  prob- 
lems of  scale  effect.  Thus,  water  tunnel  experiments  are  run  at  high 
speeds  and  the  friction  scale  effect  is  small.  The  two  remaining  require- 
ments concerning  similarity  of  advance  coefficient,  J,  and  cavitation  number. 


a,  can  easily  be  satisfied  in  the  water  tunnel  as  shown  below. 


WMIjillW 
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In  general,  the  cavitation  number  is  defined  as 


o = 


P - P 

static  vapour 

l/2pV2 


In  the  case  of  the  ship,  the  expression  above  becomes 


P - P - P 
atm  hyd  vap 


l/2pVA2 


where 


P = atmospheric  pressure 
atin 


2116  lb/ft 


Phyd  = hydrostatic  Pressure  = pgh 


P = salt  water  density 

g = gravitational  acceleration 

h = depth  to  propeller  centreline 

P = vapour  pressure  (from  tables) 

vap 


1.99  slugs/ft' 
32.2  ft/sec2 
20  ft. 

44 . 56  lb/ft2 


V = advance  speed  = (l-wT)  x V,. 


V = shiD  sneed 

s 


42.23  ft/sec 


Taylor's  wake  fraction 


0.025 


If  we  substitute  all  the  above  values  in  the  equation  for  a we  get  that 

s 

for  the  particular  ship 


I'iljllflll'ip!1'! 


’llil',l|l,J|lw’l ,"1' ' ll 
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For  the  water  tunnel  the  cavitation  number  expression  becomes 

st  vap  1 
°T  = 1/2  ptVt2  .3591 

where 

p = static  pressure  (in  mmHg) 

P = vapour  pressure  ( from  tables ) 

p = tunnel  water  density 
T 

V = tunnel  water  speed  (in  ft/sec) 

* 

Equating  the  two  cavitation  numbers  for  the  tunnel  and  the  ship,  we  obtain 
a relationship  between  and  P^  in  the  form 

Pst  - oV^  + 6 (-1) 

where  a and  B are  known  constants. 

A second  equation  can  be  obtained  by  equating  the  values  tor  J *or 
the  ship  and  the  model 

V?  x 60 

JT  = Js  = ND 


26nmHg 

1.93  sluss/ft^ 


where 

J = advance  coefficient  for  the  ship 

s 

D = model  propeller  diameter 
J = advance  coefficient  for  the  model 


.8205! 
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M = model  propeller  r.p.m. 


Thus  we  have  a second  relationship  between  the  tunnel  parameters  in  the 
form 

VT  = yN  (2) 


where  y is  again  a constant. 

A graphical  representation  of  equations  (1)  and  (2)  is  given  in 

fig-  ( 22  and  23).  Clearly,  by  choosing  a value  for  N,  values  of  V and 

P can  immediately  be  obtained  from  these  graphs.  A limitation  to  the 
st 

minimum  value  of  V is  imposed  by  the  high  Reynolds  number  requirement  and 
so  values  below  llft/sec  are  not  very  desirable.  The  final  choice  of  tunnel 
operating  parameters  used  in  the  experiment  was  based  upon  the  natural  fre- 
quency characteristics  of  the  system  as  emphasized  in  Appendix  I. 


Water  Tunnel  Wall  Interference  Correction 

Due  to  the  presence  of  the  boundary  it  is  not  possible  to  develop  free 
stream  flow  in  the  tunnel.  As  a result,  any  information  related  to  rudder 
characteristics  obtained  in  the  experiment  must  be  corrected  for  wail  inter- 
ference, prior  to  being  made  applicable  to  the  ship  itself. 

Tor  a closed  working  section  tunnel  the  condition  that  needs  to  be 
satisfied  at  the  boundary  is  that  the  normal  to  the  wall  component  of  the 
velocity  must  be  zero.  This  boundary  condition  can  be  satisfied  analytically 
by  the  introduction  of  a suitable  series  of  images  of  the  model,  the  inter- 
ference experienced  being  the  induced  velocity  corresponding  to  the  vortex 
system  of  these  images  [14]. 


j 
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Fig. 23:  Static  pressure  versus  tunnel  water  speed  to  simulate 
ship  cavitation  number. 
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In  the  case  of  a tunnel  with  rectangular  cross  section  of  breadth  ”b” 
and  height  "h"  the  upward  inclination  of  the  stream  caused  by  the  coshined 
effect  of  the  whole  system  of  images  is  given  by 


S X C mm  2 ? 2 

c = f- 1 l <-!)”  V A ; 1 s *4-  « c, 

hrtV) 


where  X - 


= total  wing  area  of  model 
= lift  coefficient 


~ cross  sectional  area  of  the  tunnel  = b * h 


-»  n 2 2 2 • 

/£(-l)  a -A  n _ *2  S*2  T 

— -i  t- 


2 2 2 2 
(m  +A  n ) 


1=1  2>*i  ' * 

1 f e 


Sufficiently  accurate  results  may  be  obtained  by  retaining  only  the  first 
tens  of  the  last  exponent  series  and,  hence. 


*2  8*2 

* " 3 + 2>» 

1 + e 


And  therefore. 


S x C 


X * = 6 X _d_  * 


h i 

« *♦  - -E-”  "S’ 


he  drag  coefficient  Cn  needs  also  to  be  corrected  for  wall 


interference 
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by  ACp  given  by 

AC  = t * C,  = 6 * - x r2  (5) 

D L C*  *• 

The  corrections  in  both  angle  of  incidence,  c,  and  drag  coefficient,  ACjj, 
are  incorporated  in  the^steady  state  data  reduction  program. (see  Appendix  VI). 


APPEKDIX  VI 


Data  Reduction  Ccrcuter  Prograns 

1.  Functional  Description  - Steady-State  Bata  Reduction  Prcgraan 


This  progress  accepts  as  input  raw  data  recorded  directly  free 


j 


580  Read  run  input  data 

NTAP  = Tunnel  velocity  nanometer  pressure  tap  identification 
NFLD  = Manometer  fluid  identification 
ANGL  = Angle  of  attack 


R = Six  component  strain  gage  dynamometer  output 
700-800  Correct  input  data  for  zero  drift 
960-1310  Non-dimensionalized  forces  and  moments 
1540-1580  Correct  non-dimensional  coefficients  for  tunnel  wall  effect 


84b 


W 


9 


| 

% 


sv 


2.  Functional  Description  - Dynamic  Data  Reduction  Program 
This  program  performs  two  functions: 

(a)  Obtains  the  dynamic  response  of  the  strain-gage  dynamometer  at 
the  test  frequency; 

(b)  Reduces  the  test  data  using  the  dynamic  calibration  obtained 
in  (a). 

A detailed  description  of  each  segment  of  the  program  is  as  follows: 

Lines  Function 

» 

180  Read  input  constants 

COIL  = Calibration  coil  constant  (lbs/ampere) 

ECAL  = Voltage  drop  across  calibrating  resistor  CRMS  volts) 
RCAL  = Calibration  resistor  value  (OHMS) 

DIAM  = Propeller  diameter 
RFM  = " RPM 

RO  = Fluid  mass  density 

Y1,Y2  = Spanwise  points  of  application  of  electromagnetic 
shaker 

KSTOP  = Control  parameter  to  stop  calculation 

200  Compute  calibrator  force  FCAL 

210  Read  response  of  dynamometer  channels  FA  and  BX  to  the 

calibrator  force  applied  at  the  two  spanwise  points  of 
application  (index  1=1  or  2).  These  are  supplied  in  the 
form  of  an  amplitude  (J=l)  and  a phase  (J=2). 


Convert  input  FZ  and  BX  amplitude/phase  representation 
to  real/ imaginary  representation  (FZMAG,BXMAG) . The 
latter  are  then  converted  to  FORTRAN  complex  variables 
CFZ,  CBX. 


250-310 
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320-380  Compute  dynamic  influence  coefficient  matrix  CA-CD  using 
built-in  FORTRAN  complex  arithmetic 

560  Force  coefficient  normalizing  factor  pn2Dl* 

620  Read  amplitude  and  phase  of  dynamometer  output  during  test 

650-700  Convert  to  complex  variables  CRUNFZ,  CRUMBX  as  in  lines 

250-310 

740-C50  Convert  input  to  forces  by  applying  dynamic  calibration 
values,  non-dimensionalize  and  convert  from  complex  to 
amplitude/ phase  representation 

920-940  Write  results  and  transfer  to  210  for  next  set  of  test 
. . data 
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00010 

00020 

00030 

00040 

00050 

00060 

00070 

00080 

00090 

00100 

00110 

00120 

00150 

00140 

00150 

00160 

00170 

00180 

00190 

00200 

00210 

00220 

00230 

00240 

00250 

00260 

00270 

00280 

00290 

00300 

00310 

00320 

00330 

00340 

00350 

00360 

00370 

00380 

00390 

00400 

00410 

00420 

00430 

00440 

00450 

00460 

00470 

00480 

00490 

00500 


C STEADY  STATE  DATA  REDUCTION  PROGRAM 
INTEGER  RS, S, BLANK 
REAL  MX, MY, MZ, LIFT, MAC 

DIMENSION  IDENTE18),ZME2),ZIE7,2),ZFE2),NTAPE 40) ,NFLOE40), ANOME 40) 
1,  ANGL(40),R( 40,  7),DZIE7),CLE40),CDE40),CME40),CPLE40),CYE40), 
2VCES),VEE5),FME2),CE7),CLDE40).CLSQE40),RNE40),SIE4),COE40) 

3, RHS  E4), ANGDE 40) 

DATA  KI/5/,KO/6/ 

C VC  AND  VE  ARE  CONSTANTS  RELATED  WITH  THE  CONVERSION  OF  THE 
C MANOMETER  READING  INTO  WATER  SPEED.  WHICH  OF  THE  5 VALUES 

C GIVEN  SHOULD  BE  USED  IS  DETERMINED  BY  THE  CHOICE  OF  THE 

C PRESSURE  TAPS  AT  WHICH  PRESSURE  WAS  MEASURED 

DATA  VC/4. 1144  7,  2.6  0088, 0.  9016,  0.6214, 4. 7381/, VE/. 49353, .5,. 4 724,0 
1.50638,  .5/ 

READ  EKI,*)TUIST, SHAFT 
104  FORMAT(6F8.5,2F12.5) 

READ  { K I , * 3 NRT, NTT, AREA, SPAN, MAC, AZL,VELINC 

100  FORMATE I4,I4,F9,2,FS.3,F9.4,F8.2,F9.4) 

I F ( AREA. LE . 0 ) GO  TO  99 

C RHO  CORRECTS  THE  VALUE  OF  THE  TUNNEL  WATER  DENSITY  FOR  CHANGES 
C IN  TUNNEL  TEMPERATURE. FM  AMD  SCALE  TAKE  CARE  OF  THE  EXPANSION  OF 

C THE  FLUID  AND  THE  MANOMETER  SCA^E  RESPECTIVELY,  DUE  TO  CHANGES  IN 

C THE  ROOM  TEMPERATURE. 

RHO -1.9574-0. 0002 8 -NTT 
SCALE-1.000952-13.3E-6-NRT 
FME1) "1 . 4875- . 00071 -NRT 
FM(2)«24.5054-. 00236-NRT 
WRITE  (KO, 200) 

200  FORMAT  E'1',//17X, 'RUDDER  INPUT  DATA’//2X, 

l'TR',2X, 'TT',4X, 'AREA', 4X, 'SPAN', SX, 'MAC ' , SX, ' AZ L* , 3X, ' VELI NC' ) 

WRI TE (KO, 150)NRT,NTT, AREA, SPAN.MAC, AZL, VELI NC 
150  FORMAT (I4,I4,F9.2,F8.3,F8.4,F8.2,F8.4) 

30  READ  (KI,103)(IDENT(N),N-1,18) 

103  FORMAT  (18A4) 

C AT  THIS  POINT  THE  MANOMETER  AND  STRAIN  GAUGES  'ZERO  READINGS', 

C <ZM  AND  Z1  RESPECTIVELY), BEFORE  AND  AFTER  A RUN  ARE  READ. IN  THE 

C SUBSEQUENT  STEPS  THE  PROGRAM  DECIDES  THE  DEGREE  TO  WHICH  THE  DIF- 

C FERENCE  BETWEEN  THE  'ZEROES’  BEFORE  AND  AFTER  A RUN,  SHOULD  AFFECT 

C THE  READINGS  TAKEN  DURING  THIS  RUN. 

REAOEKI,-)EZMEK),EZI (M, K),M-1, 6) , K-l, 2 ) 

101  FOR!!AT(F5.0,10X,5{2X,F5.0>) 

WRITE  (KO,600)(IDENT(N),N»1,18) 

600  FORM  AT ( / / , 1 X , 1 8A4 / / ) 

WRI TE ( KO, 201 ) 

201  FORMATE/,'  ZERO  READINGS  BEFORE  AND  AFTER'/’  MAN  FX  FZ 
1BX  BZ  FY  BZ’) 

WRITE(KO,  204)(ZM(K>, C Z I (M, K ) ,M-1, 6 ) , K-l, 2 ) 

204  FORMAT (2(F4.0,6(2X,F5.0)/)) 

WRITE(KO,212)TWIST,SHAFT 

212  FORMATE'  TWIST-*, F7. 1, ' SHAFT  Dl A. • ’ , F7 . 2, ’ IN.’) 
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00510 
00520 
00530 
00540  C 
00550  C 
00560  C 
00570  C 
00580 
00590 
00600 
00610 
00620 
00650 
00640 
00650 
00660 
00670 
00680 
00690 
00700 
00710 
00720 
0073C 
00740 
00750 
00760 
00770 
00780 
00790 
00800 
00810 
00820 
00830 
00840 
00850 
00860 
00870 
00880 
00890 
00900 
00910 
00920 
00930 
00940 
00950 
00960 
00970 
00980 
00990 
01000 


00  2 J-l, 40 

JT-J-1 

JJ-J 

NTAP  AND  NFLD  ARE  THE  CODE  NUMBERS  OF  THE  TAP  AND  FLUID  CHOSEN 
AHOM(J)  IS  THE  MANOMETER  READING  AT  SOME  ANGLE  OF  ATTACK  DENOTED 
BY  THE  INDEX  J.  R!J.M)  ARE  THE  OUTPUTS  OF  THE  STRAIN  GAUGES  AT 
THE  SAME  ANGLE. 

READOCI  , •) NTAP! J), NFLD! J), ANOM!  J), ANGL! J) , !R!J,M), M-1,6) 

102  FORMAT! 1 4, 1 1,F5 . 0,  F7.  2,  2X, F8 . 2, F8 . 2, F8 . 2, F8 . 2, F8.  2,  F8. 2 ) 
IF!ANOMiJ).LE.0.0)  GO  TO  3 
2 CONTINUE 

203  FORMAT! IX, 1 2, 1 1,F7. 0,  F5 . 1, 2X, F6. 0, 3X, F6.0, 3X, F6.0, 3X, F6.0, 3X,F6. 0, 
13X,F6.0) 

GO  TO  99 

3 WRIT£!KO,202) 

202  FORMAT!  /'  INPUT  DATA  AS  RECORDED'/ ' TF  AHUM  ANGLE  FX 
1 FZ  BX  BZ  FY  BY  ') 

WRITE!KD,203)!NTAP!J),NFLD!J),ANOM!JJ,ANGL!J),!R..  J,M), M-1,6),  J-l,  J 
1J) 

BUG-1. 0/!JT-l) 

DZM-! ZM!2)-ZM! 1 ) )-RUG 
DO  4 M-1,6 

4 DZI !M)-!ZI !M, 2)-ZI !M,  1) )«BUC 
DO  7 J-1,JT 

IF! J.EQ.l)  GO  TO  6 
IFtNTAP(J).EO.O)  NTAP! J) -NTAP! J-l) 

IF!NFLD!J).EQ.0)  NFLD! J) -NFLD! J-l) 

6 BUG-J-1 

ANOM! J ) -ANOM! J)-ZM! 1 ) -BUG-OZM 
ANGL! J) -ANGL! J)-AZL 
DO  7 M-1,6 

R!J,M)-R!J,M)-ZI  !M,  D-BUG-DZI  (M) 

7 CONTINUE 

WRI TE(KO, 205)  ( I OENT!N), N-l,  18 ) 

205  FORMAT  (.lilt  IX,  18A4/  / ) 

WRI TE !KO, 705) !NTAP( J), NFLD! J), ANOM! J), ANGL! J),(R! J,M), M-1,6 ), J-l, J 
IT) 

705  FORMAT!'  INPUT  DATA  CORRECTED  FOR  ZERO  READINGS  '/!I3, 

1I1,F6.0,F7.2,2X,6F9.2)) 

33  WRITE  (rO,207)!IOENT!N),N-l,18) 

207  FORMAT!///, IX, 18A4//) 

607  WRITE! KO, 707) 

707  FORMAT t ' RUDDER  DATA  REDUCTION  ••'//'  ALPHA  FX-LB 

1 FZ-LS  BX-INLB  8Z-INLB  ORAG-LB  LIFT-LB  VER-LB  MX- INLB  M 
2Z-INLB  MY-INLB  VEL-FPS') 

BEAR-MAC/ ! 12*3. 9739*tXP!67. 6832/NTT)) 

DO  9 J-l, JT 
I -NTAP! J) 

IFO.LT.l.OR.I.GT.S)  GO  TO  99 
K-NFLDl J) 


s 
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01010  IF(K.LT.1.0R.K.GT.2)  GO  TO  99 

01020  BUG*VE(I ) 

01030  C THE  WATER  SPEED  IS  COMPUTED  AT  THIS  POINT  FROM  THE  KNOWLEDGE 

01040  C OF  THE  VARIOUS  CONSTANTS  RELATED  TO  THE  MANOMETER  FLUID  AND 

01050  C THE  FACTOR  CORRECTING  FOR  THE  INTRODUCTION  OF  THE  SPLITTER 

01060  C PLATE  IN  THE  FLOW  (VELINC). 

01070  V*(ANOM(J)*FM(K)*SCALE/(VC(l )*RHO)  )**8UG»(  1. 0*VEL1  NC> 

01080  C AT  THIS  POINT  THE  PROGRAM  MAKES  USE  OF  THE  STATIC  CALIBRATION 

0109C  C RESULTS  TO  CONVERT  THE  FORCES  I Ni  l LBS.  AND  THE  MOMENTS  INTO 

01100  C IN-LBS  TAKING  ALSO  CARE  OF  THE  OOSS-TALK  BETWEEN  THE  DIFFERENT 

OHIO  C GAUGES. 

01120  FX— 0.076768795*R(J,l)-0.06336i>/ll*R<J,4) 

01130  FZ "0.07236 94 65 *R(J,t)  + 0. 137756 *R(J,  3) 

01140  BX »- 0. 051182352 »R(J# 2 )*0. 4552055 23 »R(J, 3) 

01150  BZ «-0. 042513862 •R(J,1)*0. 22310536 »R(J, 4) 

01160  FY-0.0 

01170  BT—  0.6493S*R(J,6) 

01180  RAT-(ANGL<J)*AZL) *0.0174532 

01190  SKJ)-SIN(RAT) 

01200  CO(J)«COS<RAT) 

01210  C HERE  THE  LOADS  MEASURED  ARE  RESOLVED  IN  THE  LIFT  AND  DRAG  DIREC- 
01220  C TIONS.  THE  TUNNEL  WALL  INTERFERENCE  CORRECTION  HAS  NOT  YET  BEEN 
01230  C MADE. 

01240  LIFT -FZ*CO( J)-FX*SI ( J) 

01250  DRAG*FX*CO<J)*FZ*SI(J) 

01260  MX-BX*CO(J)*BZ*SI(J) 

01270  MZ-BZ*CO<J)-8X*SI(J) 

01280  VER-FY 

01290  MY-BY 

01300  ANGO ( J)  -ANGL ( J )*MY/TWI  ST 

01310  WR|TE<KO,206)ANG0(J),FX,FZ,BX,BZ,DRAG,UFT,VER,MX,MZ,MY,V 

01320  206  F0RMAT(11F9. 2, F8. 2) 

01330  R VS *0. 5*RHO • AREA* V**2/ 144.0 

01340  C HERE  THE  DIFFERENT  FORCES  AND  MOMENTS  ARE  NON-DIMENSIONAL  I ZED 

01350  CL(J)-LIFT/*VS 

01360  CLSQ(J)*CL(J)**2 

01370  CD(J)-DRAG/RVS 

01380  CLD(J)-CL(J)/CD(J) 

01390  CM( J ) »MY/R VS/MAC 

01400  CPL(J)-(MX/LIFT)/SPAN 

01410  CY(J)-VER/RVS 

01420  9 RN(v)»V*BEAR 

01430  WRITE  (KO, 208) ( I DENT(N), N»l, 18) 

01440  208  FORMATC// / , IX, 13A4/ / ) 

01450  608  WRI TE(KC, 708) 

01460  708  FORMAT(*  «*RUDDER  DATA  IN  NON-DIMENSIONAL  FORM**'//'  ALPHA 

01470  2CL  CD  C"  CPL  CY  L/D  RN*10**-6  CLSO’ ) 

014P  WR I T£(KO, 209) ( ANGD( J),CL( J), CD( J),CM{ J),CPL( J),CY( J), CLD( J),RN( J), 

01490  1CLSQ<J),J*1,JT) 

01500  209  FORMAT (F8 .2,F7.3,2F8.4,3F7.3,F8.3,5X,F7.4) 
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01510  00  15  J-l.JT 

01520  C HERE  THE  RUDDER  ANGLE  AND  THE  DRAG  COEFFICIENT  ARE  BOTH 

01530  C CORRECTED  FOR  TUNNEL  WALL  IN’ERFERENCE. 

01540  ANGD(J)»ANGD{J)*0. 02045 *AREA»CL(J) 

01550  CD( J)*CD( J)*0 . 000S570*AREA»CL< J)*»2 

01580  15  CLD(J)-CL(J)/CD(J) 

01570  WRI  TEOCO, 211)  < IDENT(N),?'«1. 18) 

01580  211  FORMATC/ //, 1X# 18A4// ) 

01600  711  FORMATS 'nl>PRIOP  DATA  CORRECTED  FOR  TUNNEL  INTERFERENCE'//*  AL 

01610  1PHA  CL  CD  CD  CPL  CY  L/D  RN*10«*-6  CLSO 

01650  2 WRI  TEUCO,  209  ){ANGD<J),CL{J),CD<J),CM(  J),CPL(J),CY(J),CLD(J),RN(J), 

01640  1CLSQ(J),0-1,JT) 

01650  GO  TO  30 

01660  99  STOP 

01670  END 
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RUDDER  ItIPUT  DATA 

TR  TT  AREA  SPAN  MAC  AZl  VELINC 

82  85  74.19  9.354  7. 9670  0.0  0.0910 


THE  FOLLOWING  RESULTS  ARE  FOP  A CLEARANCE  RATIO  OF  .37 
ZERO  READINGS  BEFORE  AND  AFTER 


:an 

FX 

FZ 

8X 

BZ  FY  BZ 

0. 

0. 

500. 

-650. 

0.  0.  0. 

0. 

-56. 

494. 

-648. 

-4.  0.  0. 

TWIST* 

10000.0 

SHAFT  DIA. * 2.00  IN. 

l 

TF 

INPUT  DATA  AS 
AHOIt  ANGLE 

RECORDED 

FX 

FZ 

BX 

BZ 

FY 

BY 

11 

1506.  0.C 

-16. 

622. 

-666. 

-f  8 . 

0. 

-3? 

00 

1506.  5.0 

-19. 

90. 

-24. 

-67. 

0. 

-70 

00 

1506.  10.0 

-65. 

-440. 

632. 

36. 

0. 

-58 

00 

150.'.  15.0 

-154. 

-97?. 

1318. 

258. 

0. 

-102 

00 

1506.  26. 0 

-286. 

-1504. 

1964. 

590. 

0. 

-96 

00 

971.  25.0 

-278. 

-1158. 

1475. 

67C. 

C. 

-52 

00 

S71.  30.0 

-444. 

-1472. 

1256. 

1020. 

0. 

-29 

00 

-100.  0.0 

0. 

0. 

0. 

C. 

0. 

0 

THE  FOLLOWING  RESULTS  ARE  FOR  A CLEARANCE  RAT|D  OF  .37 
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1506. 

0.0 
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122.00 

-16.00 

-C8.0G 

0.9 

-38.00 

11 

1506. 

5.00 

-13.00 

-409.00 

625.67 

-66.33 

0.9 

-70.00 

11 

1506. 

10.00 

-53.00 

-938.00 

1281.33 

37.33 

0.-'' 

-38. 0C 

11 
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15.00 
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-1475.90 

1967.09 

2C0.C" 

O.o 

-102.9" 

11 

1506. 

20.00 

-262.00 

-2000.00 

2C12.C7 

5P2.67 

0.9 

-96.D5 

11 

971. 

25.00 

-248.00 

-1F53.0O 

2!" 

'7  .33 

-52.99 

11 

971. 

30.00 

-408.00 

-1956.0'' 

2504.00 

1924.00 

o!o 

-20. CO 
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THE  FOLLOWING  RESULTS  ARE  FOR  A CLEARANCE  RATIO  OF. 37 
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ZERO  READINGS  BEFORE  AND  AFTER 

MAN  FX  FZ  BX  BZ  FY  BZ 
0.  0.  -S06.  530.  0.  0.  0. 

0.  -2$.  -4f2.  630.  A.  0.  0. 

TWIST-10000.0  SHAFT  Dl A. - 2.00  in. 


INPUT  DATA  AS  RECORDED 


TF 

AIIOM 

ANGLE 

FX 

FZ 

BX 

BZ 

FV 

BY 

11 

1506. 

0.0 

-3. 

-A00. 

616. 

-72. 

0. 

-37. 

00 

1506. 

-5.0 

-26. 

126. 

6. 

1*. 

0. 

-r. 

00 

1506. 

-10.0 

-BA. 

698. 

-f  34. 

204. 

0. 

12. 

00 

1506. 

-15.0 

-156. 

1312. 

-1320. 

464. 

0. 

12. 

00 

1506. 
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-210. 
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-2020. 

714. 

0. 

-2. 

00 

971. 

-25.0 

-22' . 

1432. 

-1402. 

724. 

0. 

-24. 

00 

971. 

-30.0 

-23A. 

2035. 

-2032. 

870. 

0. 

-r,o. 

00 

-100. 

0.0 

0. 

0. 

0. 

0. 

0. 

0. 

THE  FOLLOWING  RESULTS  AP.E  FOR  A CLE77A  XE  OF.37 


INPUT  nATA  CORRECTED  FOR  ZERO  REA^IJr- 


11 

1506. 

0.0 

-3.00 

106.00 

-14.00 

11 

15SC. 

-5. CO 

-21.33 

628.00 

-624.30 

11 

1506. 

-10.00 

-74. F7 

me. on  -1264.00 

11 

1500. 

-15.30 

-142.00 

1806.00  -1050. on 

11 

1506. 

-20.00 

-191.33 

2488.00  -2650.00 

11 

371. 

-25.00 

-202.67 

1978.00  -2122.00 

11 

971. 

-30.00 

-206.00 

2517.00  -2662. on 

i 


-72.DC 

0. 

-37. C* 

17.33 

o.  0 

-6.?C 

202.67 

»■•!? 

12.  or 

4G2.O0 

O.o 

12. CO 

711.33 

o.O 

-2.C3 

720.67 

p..o 

-24.03 

866.00 

O.o 

-60.00 

Table  III  (con’t) 
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DYNAMIC  DATA  REDUCTION  PSOG^A" 

THE  FIRST  PART  OF  THIS  PROGRAM  IS  USED  TO  CALIBRATE  THE 

TWO  DYNAMOMETER  CHANNELS  HE  A s£  EXAMINING. A (V3M  FORCE  *FCAL?, 

IS  APPLIED  ARP  THE  RESPONSE  OF  EACH  CHARNEL  TO  THIS  FORCE, *FZ*. 

A*  ’SI*  RESPECTIVELY,  IS  COMPUTED  AT  TmO  CIFFEPENT  LOCATIONS 
OF  THE  APPLIED  FORCE, THESE  RESULTS  CAR  RE  CONFINED  TO  GIVE  AN 
INFLUENCE  COEFFICIENT  MATRIX  CONTAINING  THE  ELErENTS  •CA’.'C**, 

* CC*, * CP', WHICH  ARE  USED  LATER  ON  DURING  THE  EXPERIMENT  TO 
EVALUATE  THE  AMPLITUDE  AND  PHASE  ANGLE  OF  THE  TVO  RESPONSES, 

AS  INDICATED  IN  THE  SECOND  PAST  CF  THE  PROGS  Aft. 
complex  CFZ(j),CRi{2>,MUNFz.cauNtx,CA,c9,ce,eo,er9RCE,CHPix,o*arN 
IT.CDENOM 

DIMENSION  FZC2,  21.RXC2,  21,FZ5tASC2,  ?),9XMA0C 2,  21  »RUWNXC2),IRJICFZC  21 
E,RFZC2},R9X(2],LA3ElCltl 
Data  xi/s/,xo/e/ 

HR|T£(S,£0) 

FORMATE!!,  'EXECUTION  REGINS’//) 

READ  CXI,-lCOIL,ECAl,RCAL,OIA!t.RPW,R3,YI,V2.XSTOP 
IF  CKSTOP.EQ.Cl  GO  TO  S3 
FCAL*COi l-ECAl-l.LlL2Ik/RCRL 

REAPCXl  ,•)  {CFZ£t,Jl.J-1.21.l-I,2>,CC9ICl,J>,J-;,21,l-I»21 
NEXT  THE  PROGRAM  READS  THE  CALIBRATION  VALUES.  THE  FIRST  INDEX 
REFERS  TO  LOCATION  ON  THE  RLPOER  AND  THE  SECOND  IF  1 IS  THE  AMPLI- 
TUDE OF  THE  SIGNAL  AND  IF  2.1$  THE  PHASE. 

03  IS  X-1,2 

FZMAS{K,l)-FZ<K,l)*CSS(B.CmSJ2B*FZ(K.2>) 

»XWRSCK,X)-$XCK,l!*CGSCI.S17fcS329»9XEK,2>} 

FZHA5EX , 21 -FZCK, 1 1-SINC3. S12AS JJS-FZf R» 2) > 
SK»AG(K.2)-«X(K,l)*SIN(6.«i7t5m-SX{K,2]> 
CFZ(K}-CMPLI(FZMAG(K,1>,FZMAT.(K,2I} 

CSX  C X 1 -CHPII  (RIM  AG  EX,  J 1 , R IMAGE  R,  2 1 1 
C5EWW.csiCll-CF2t21-CBXt21-CFZEIl 

AT  THIS  POINT  TW  PROGRAM  HAS  REWRITTEN  THE  INPUT  R FACINGS  IN 
COMPLEX  FORM  AMD  IS  READY  TO  COMPETE  THE  INFLUENCE  COEFFICIENTS 
CA-((CBI(1)-CBX(2})/CSENPM)*FCAL 
CB-({CFZf2l-CFZtll 1/CDENflMl-FCAl 
CC-C  CCIXC 1 1 -T2-CRX { 21 *Y» J/C3EN0MJ-FCAL 
CD-C  CCFZ  C2 1 -Yl-CFZf 1 J -Y2 1/ CPC WOW  » -FCAL 
WRITE  (XD  2R1  C0ll,ECAl,RCAl,Dl*“,RPM,R0,YI.Y2,FCAl 
FORMATC'l’.M.’Ceil-’.FS.J.Xl.’ClS/AMPl'./.Sl.'CCAl-'.rS.J.X 
XX,  * C VOLTS 1*,/,5X, *RCAl— *,FS  . J,  IX,  “COIR'S I * ,/,SX,  "PROP. SIAM** 

2, FI. S, IX,  * C INI  *,/,SX,  *RBi«*,F|.  l,/,5*,  ,DENSI7T»,,F1R.S, If, 

S’ C L9S/FT--JJ  *,/, Si, *LI-’,FS. S,IX, *f INJ',/,51, *L2-*,F$. J, IX, 
*,CIN)*,///,SS.,CAllSRATION  FORCE-*. Fie.S,lX»'CUSl*l 
FORMATE///. 2RI.  'CALIBRATION  REAPINGS*,///,  2SX,  *A»tPl*,fX,  ’PHASE*, 
XIJX.'AMPt*,*X, 'PHASE*. //.2X. 'FORCE  AT  LX ', 2X, *FZR£AP-’,F 
2 13. S, LX, FI .1,21, 'IXREAD-'.FXR.S, LX.FB.L, /, 2X, ’FORCE  AT  L2*. 
J2X,*FZREAO-*,F10.S.LX,F».L,2X,**XREAP-*,F13.S,LI,FI.L/J 
WRITECX0.2I1  EFZEl,ll,FZ(l,21,Rltl,ll,RXEl,21, 1-1,21 
WRITEEX0.221  EFZMA&EI,X},FZMAGEI,2},BXMAGEI,I},>XHAGEI,21 


;arle  .V: 


jdder  dvnacic  data  reduct 


ion  program 


00510 

00520 

00530 

00540 

00550 

00560 

00570 

00580 

00590 

00600 

00610 

00620 

00630 

00640 

00650 

00660 

00670 

00680 

00690 

00700 

00710 

00720 

00730 

00740 

00750 

00760 

00770 

00780 

00790 

00800 

00810 

00820 

00830 

00840 

00850 

00860 

00870 

00880 

00890 

00900 

00910 

00920 

00530 

00940 

00950 

00960 

00970 

00980 


ii  f mu 
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1. 1-1.2) 

22  FORMAT*///, IX, 'AT  L1',3X, ' REALFZ-* , F10. 5, IX, ' I MAGFZ* ' , F10. 5, 

11X, 'REALBX-',F10.5,1X, ' I MAGBX- ' , F10. 5, //, IX, 'AT  L2',3X, 
2'RCALFZ«',F10.5,1X, ' IMAGFZ-', F10. 5, IX, ' REALBX-', F10.5, IX, 

3 ' IMAGBX« ', F10. 5/ ) 

BUG-12. 0**4/ (RO*( (RPM/ 60. 0)»*2)»(DIAM«*4)) 

2 READ* Kl ,40)  ( LABEL (N) , N*l, 18 ) 

40  FORMAT (18A4 ) 

C AT  THIS  POINT  THE  PROGR\M  READS  THE  VALUES  OF  THE  LOADS  MEASURED 

C IN  THE  ACTUAL  EXPERIMENT.  THE  INDEX  AS  BEFORE  DIFFERENTIATES 

C BETWEEN  AMPLITUDE  AND  PHASE. 

REAPOCI,*)  RFZ(1),RFZ(2),RBX(1),RBX(2), NCODE 
42  FORMAT (4F10.5, 110) 

IF  (NCODE. EQ.O ) GO  TO  1 

"UNBX (1 )«RBX(1 ) *COS(0. 01745329*RBX(2) ) 

RUNBX(2)»hBX(1)*SIN(0.01745329*RBX(2)) 

RUNFZd ) ■RFZ(1)*COS(0.01745329*RFZ(2)) 
RUNFZ(2)*RFZ(1)*SIN(0.01745329«RFZ(2>) 

C RUNFZ«Ct1PLX(kUNFZ(l),  RUNFZd) ) 

CRUNBX-CMPLX(R'JNBX(1),RUNBX*2)) 

C CFORCE  AND  CMOMNT  ARE  THE  FORCE  AND  MOMENT  VALUES  OF  THE  READINGS 

C CORRECTED  FOR  THE  CROSS  TALK  BETWEEN  THE  TWO  GAUGES.  THEY  ARE  IN 

C COMPLEX  FORM  AND  THEIR  AMPLITUDES  ARE  IN  LBS. 
CFORCE-CRUNFZ-CA+CRUNBX-CB 
CMOMNT“CRUNFZ*CC*CRUNBX*CD 
FREAL "REAL (CFORCE) 

F I MAG-AI MAG (CFORCE) 

BREAL "REAL (CMOMNT) 

BlMAG-AIMAG(CMOMNT) 

C FAMP  AND  BAMP  ARE  THE  FORCE  AND  MOMENT  AMPLITUDES  IN  NON-DIMENSIO- 

C NAL  FORM.  THEY  CORRESPOND  TO  'KF-1000'  AND  'KB-1000'  ON  THE  GRAPHS 

FAMP-SQRT(FREAL*-2+F IMAG--2 )*BUG*1Q00. 0 
BAMP“SQRT(BREAL**2+BIMAG**2 )*BUG»1000. 0/DI AM 
FPHASE-ATAN(F IMAG/FREAL)*5  7. 2958 
BPHASE»ATAN(B IMAG/BREAL) »57. 2958 

46  FORMAT*//, IX,  18A4,//,  28X, ' EXPERIMENT  READI  NGS’ , III . 20X, ' AM. PL ' , 

16X, 'PHASE',  23X,  ' AMPL ',  6X,  ' PHASE ',//,  5X,  ’RFZ',7X,F10.5,3X, 

2F8 . 4, 7X, 'RBX',7X,F10.5,3X,F8.4,//,5X, 'FREAL- ' , F10. 5, //, 

35X, 'FIMAG-',F10.5,//,5X, ' BREAL- ' , F10. 5, //, 5X, ' Bl MAG- ' , F10. 5, / 
4/,5X, ' F AMP-1000*', F10. 5, 9X, ' FPHASE-' , F10, 5, //, 5X, 'BAMP-1000-' 
5,F10.5,9X, ' B PHASE «', F10 . 5 ) 

WRI TE(KO, 46)(LABEL(N),N»1,18),RFZ(1),RFZ(2),RBX(1),RBX(2), 

1FREAL, F IMAG, BREAL, B I MAG, FAMP, F PHAS  E, BAMP, BPHASE 
GO  TO  2 

50  WRI TE(KO,  55) 

55  FORMAT  Cl') 

STOP 

END 


Table  IV  'oon't) 
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COtL-O. 750  (LB/AMP) 
ECAI.-4.000  (VOLTS) 
RCAL-4.600  (OHMS) 
PROP.OIAM-  9.846  (IN) 
RPM-1500.0 
DENSITY- 
Ll-3 . 000  (IN) 

L2-9.000  (IN) 


1.93780  ( SLUOS/FT -»3) 


CALIBRATION  FORCE* 


0.92231  (LBS) 


CALIBRATION  READINGS 


A’lPL 

PHASE 

AM  PL 

PHASE 

FORCE 

AT 

LI 

FZREAD- 

0.2C56J 

117.0000 

BXREAD- 

0.12500 

270.0000 

FORCE 

AT 

L2 

FZREAD* 

0.U6O00 

90.0000 

BXREAD- 

0.60000 

270.0000 

AT 

LI 

REALFZ- 

-0.09307 

IMAGFZ* 

0.18266 

REALBX* 

-0.00000 

IMAGBX- 

-0.12500 

AT 

L2 

REALFZ- 

0.00000 

IMAGFZ- 

0.06000 

REALBX- 

-0.00000 

IMAGBX- 

-0.60000 

'lable  V:  Sample  output  of  dynamic  data  reduction  program 


t 


RUDDEP  AT  1 INCH  OFF  CENTRELINE  (DESIGN  POSITION).  CLEARANCE  RATIO  .33 

EXPERIMENT  READINGS 


AMPL 

PHASE 

AMPL 

PHASE 

RFZ  0.32000 

284.0000 

RRX 

0.65000 

36.0000 

FREAL*  -1.72153 

FIMAG*  0.87087 

RREAL*  -6.8S750 

BIMAG*  7.36398 

FAMP»1000=  3.51471 

FPHASE* 

-26.83350 

BAMF*1000*  1.86689 

B PHASE* 

-46.87341 

RUDDER  AT  1 INCH  OFF  CENTRELINE  (DESIGN  POSITION).  CLEARANCE  RATIO  .37 

EXPERIMENT  READINGS 


AMPL 

PHASE 

AMPL 

PHASE 

RFZ 

0.38500 

252.0000 

RBX 

1.15000 

45.0000 

FREAL* 

-2.12280 

FIMAG* 

1.88896 

RREAL- 

-12.49031 

BIMAG* 

12.15800 

FAMP»1000*  5.17670 

FPHASE* 

-41.66414 

BAMP*1000*  3.22514 

BPHASE* 

-44.22762 

♦ 


i 

f 


Table  V (con’t) 
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RUDDER  AT  1 INCH  OFF  CENTRELINE  (DESIGN  POSITION).  CLEARANCE  RATIO  .41 

EXPERIMENT  READINGS 


AMPL 

PHASE 

AMPL 

PHASE 

RFZ  0.12500 

288.0000 

RBX 

0.24000 

180.0000 

FREAL-  -0.41909 

FI MAG*  -0.20247 

BREAL-  -0.59555 

B IMAG*  -3.08578 

F AMP-1000-  0.84793 

FPHASE* 

25.78540 

BAMP-1000-  0.58113 

BPHASE- 

79.06931 

RUDDER  AT  1 INCH  OFF  CENTRELINE  (DESIGN  POSITION).  CLEARANCE  RATIO  .45 

EXPERIMENT  READINGS 


AMPL 

PHASE 

AMPL 

PHASE 

RFZ  0.31500 

252.0000 

RBX 

0.95000 

360.0000 

FREAL-  -0.98530 

FIMAO-  2.00929 

BREAL-  -1.40015 

BIMAG-  13.92307 

F AMP-1000-  4.07690 

FPHASE- 

-63.87798 

BAMP-1000-  2.58915 

BPHASE- 

-84.25746 

Table  V (con't) 
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RUDDER  AT  1 INCH  OFF  CEMTPELINR  (DESIGN  POSITION).  CLEARANCE  RATIO  .49 

EXPERIMENT  READINGS 


AM  PL 

PHASE 

AM  PL 

PHASE 

RFZ 

0.34500 

243.0000 

RBX 

0.90000 

18.0000 

FREAL* 

-1.22266 

FI  MAG* 

2.05548 

BREAL* 

-4.97760 

BIMAG* 

12.89301 

FAMP-1000*  4.35702  F PHASE*  -59.25436 

BAMP-1000*  2.55717  BPHASE-  -68.88991 


RUDDER  AT  1 INCH  OFF  CENTRELINE  (DESIGN  POSITION).  CLEARANCE  RATIO  .53 

EXPERIMENT  READINGS 


AMPL 

PHASE 

AMPL 

PHASE 

RFZ  0.45500 

243.0000 

RBX 

1.32500 

18.0000 

FREAL-  -1.68770 

FIMAG*  2.86159 

PRCAL-  -7.16850 

BIMAG-  18.74754 

F AMP-1000*  6.05233 

F PHASE* 

-59.46878 

BAMP-1000-  3.71375 

BPHASE* 

-69.07466 

Table  V (con't) 
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RUDDER  AT  1 INCH  OFF  CENTRELINE  (DESIGN  POSITION).  CLEARANCE  RATIO  .57 

EXPERIMENT  READINGS 


AMPL 

PHASE 

AMPL 

PHASE 

RFZ 

0.52500 

252.0000 

RBX 

1.70000 

54.0000 

FREAL- 

-3.19651 

FIMAG* 

2.39540 

8REAL* 

-20.56535 

BIMAG* 

15.04536 

FAMP*1000*  7.27701 

FPHASE* 

-36.84721 

BAMP*1000“  4 . 71474 

BPHASE* 

36.18871 

RUDDER  AT  1 l»CH  OFF  CENTRELINE  (DESIGN  POSITION).  CLEARANCE  RATIO  .61 

EXPERIMENT  READINGS 


AM«»L 

PHASE 

AMPL 

PHASE 

RFZ  0.39000 

270.0000 

RBX 

1.3500C 

72.0000 

FREAL-  -2.91310 

FIMAG-  0.97610 

BREAL*  -19.09775 

BIMAG*  6.24727 

FAMP»1000*  5.59702 

FPHASE* 

-18.52466 

BAHP*1000-  3.71787 

BPHASE- 

-18.11400 

Table  V (con't) 
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RUDDER  AT  1 

INCH  OFF  CENTRELINE  (DESIGN 

POSITION). 

CLEARANCE  RATIO 

.65 

experiment  readings 

AM  PL 

PHASE 

AMPL 

PHASE 

RFZ 

D. 19500 

261.0000 

RBX 

0.62500 

99.0000 

FREAL* 

-1.30840 

f ihag* 

0.21919 

BREAL* 

-9.03109 

B 1 MAG* 

-0.82759 

FAMP*1000*  2.41684 

FPHASE 

* -9.51000 

BAMP«1000*  1.68164 

BPHASE 

- 5.22434 

RUDDER  AT  1 INCH  OFF  CENTRELINE  (DESIGN 

POSITION). 

CLEARANCE  RATIO  .69 

EXPERIMENT  READINGS 

AMPL 

PHASE 

A»'PL  PHASE 

RFZ  0.12500 

252.0000 

RRy 

0.23500  27.0000 

FRFAL*  -0.48947 

FI  MAG*  0.57579 

BREAL*  -1.93851 

BIMAG*  3.25764 

FAMP*1000*  1.37676 

FPHASE 

• -4Q. 63249 

BAM^IOOO*  0.70140 

BPHASE 

• -59.24458 

Table  V (con't) 
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RUDDER  AT  1 INCH  OFF  CENTRELINE  (DESIGN  POSITION).  CLEARANCE  R ,TI0  .73 

EXPERIMENT  READINGS 


AMPL 

PHASE 

AMPL 

PHASE 

RFZ  0.23000 

225.0000 

RBX 

0.55000 

18.0000 

FREAL-  -0.54299 

FIMAG-  1.43159 

RREAL-  -2.75169 

BIMAG=  8.12841 

FAMP-1000-  2.78935 

FPHASE- 

-69.22874 

BAMP-1000-  1.58782 

B PHASE- 

-71.29758 

RUDDER  AT  1 INCH  OFF  CNTRELINE  (DESIGN  POSITION).  CLEARANCE  RATIO  .77 

EXPERIMENT  READINGS 


A*'PL 

PHASE 

AMPL 

PHASE 

RFZ  0.30000 

243.0000 

RBX 

0.73800 

45.0000 

FREAL-  -1.360R3 

FIMAG-  1.46830 

BREAL-  -S. 01360 

BIMAO-  8.16617 

FAMP.1000-  3.64723 

FPHASE* 

-47. 17331 

BAMP-1000*  2.11696 

BPHASE* 

-45.54025 

'able  V (con’t) 
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RUDDER  AT  1 I f'CH  OFF  CENTRELINE  (DESIGN  POSITION).  CLEARANCE  RATIO  .81 

EXPERIMENT  READINGS 


AMPL 

PHASE 

AMPL 

PHASE 

RFZ 

0.27000 

261.0000 

RBX 

0.76400 

81.000 

FREAL* 

-1.73683 

FI MAG* 

0.62921 

. 

BREAL* 

-11.25936 

BIMAG* 

2.28652 

FAMP*1000*  3.36536  F PHASE*  -19.91405 

BAMP«1000*  2.12582  B PHASE*  -11.47938 


RUDDER  AT  1 INCH  OFF  CENTRELINE  (DESIGN  POSITION).  CLEARANCE  RATIO  .85 

EXPERIMENT  READINGS 


AMPL 

PHASE 

AMPL 

PHASE 

RFZ 

0.18000 

270.0000 

RBX 

0.60000 

95.0000 

FREAL* 

-1.30647 

FIMAG* 

0.15409 

BREAL* 

-8.82014 

BIMAG* 

-0.39026 

FAHP«1000=  2.39661  FPHASE*  -6.72673 

BAMP*1000*  1.63357  RPHASE*  2.53347 


Table  V (con't) 
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RUDDER  AT  1 INCH  OFF  CENTRELINE  (DESIGN  POSITION).  CLEARANCE  RATIO  .93 

EXPERIMENT  READINGS 


AMPL 

PHASE 

AM  PL 

PHASE 

RFZ 

0.12000 

252.0000 

RBX 

0.43800 

117.0000 

FREAL* 

-0.74505 

FIMAG- 

0.01817 

BREAL* 

-5.60540 

R IMAG« 

-2.29081 

FAMP*1000*  1.35772 

F PHASE* 

-1.39727 

RAMP»1000“  1.12042 

BPHASE* 

22.22878 

RUDDER  AT  1 INCH  OFF  CENTRELINE  (DESIGN  POSITION).  CLEAR/ NCE  RATIO  1.01 
l EXPERIMENT  REAMINGS 

i 

Tt 


• 

AMPL 

PHASE 

AMPL 

PHASE 

RFZ  0.20000 

243.0000 

RBX1 

0.49000 

81.0000 

FREAL*  -1.02251 

1 

i 

FI MAG*  0.62692 

( 

RREAL*  -7.09139 

i 

BIMAG*  1.78392 

; 

FAMP-1000*  2.18504 

FPHASE* 

-31.51349 

- 

BAMP*1000*  1.35298 

BPHASE- 

-14.12039 

Table  V (con’t) 


RUDDER  AT  1 INCH  OFF  CENTRELINE  (DESIGN  POSITION).  CLEARANCE  RATIO  1.09 

EXPERIMENT  READINGS 


ANPL 

PHASE 

AH  PL 

RFZ  O.23S0O 

315.0000 

RRX 

o.socoe 

FREAL-  -1.37119 
FI HAG-  -0.76190 
BREAL-  -7.60370 
BIMAG-  -5.19159 
FA.*  >* 1000 « 2.I57S1 

FPHASE- 

29.06151 

BAMP-1000-  1.70359 

B PHASE- 

39.32390 

PHASE 

126.0000 


Table  V (con't) 


